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Abstract. We compute averages of products and ratios of characteristic polynomials associ- 
ated with Orthogonal, Unitary, and Symplectic Ensembles of Random Matrix Theory. The pfaf- 
fian/determinantal formulas for these averages are obtained, and the bulk scaling asymptotic limits 
are found for ensembles with Gaussian weights. Classical results for the correlation functions of 
the random matrix ensembles and their bulk scaling limits are deduced from these formulas by a 
simple computation. 

We employ a discrete approximation method: the problem is solved for discrete analogues 
of random matrix ensembles originating from representation theory, and then a limit transition 
is performed. Exact pfaffian/determinantal formulas for the discrete averages are proved using 
standard tools of linear algebra; no application of orthogonal or skew-orthogonal polynomials is 
needed. 
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1. Introduction 

> 

| 1.1. The problem. Consider the linear space of real symmetric, Hermitian, or quaternion real 

Hermitian square matrices H with the Gaussian measure 

^ ; (1.1.1) P(dH) = const • exp(-Tr# 2 ) dH 

These probability spaces are the basic objects of interest of Random Matrix Theory (RMT, for 
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short); they are known as Gaussian Orthogonal, Unitary and Symplectic ensembles, respectively 
(GOE, GUE, and GSE). The goal of this paper is to study the averages of products and ratios of 
characteristic polynomials 



s / det(ai - H) ■ ■ • det(a fc - H) 

l "' J \det(^-H)---det(p m -H) 

with respect to these ensembles and their generalizations. 

Despite the fact that the Gaussian ensembles have been extensively studied, the progress on 
evaluating averages (|1. 1.2(1 remained rather limited until very recently. Starting from 1995 there 
appeared a number of papers by different authors where the case of unitary ensembles (exact 
definitions are below) was essentially settled both for the random matrices of finite size and for 
their scaling limits in the bulk of spectrum as the size of matrices tends to infinity, see Refs. 
IH US! US! SSI EH E3 ESI EOl ESI ISH [HI HI [HI EZl The results turned out to be closely related to 
some problems of the classical number theory, see Refs. |2l)j . [SH1-|121 f° r details. 

Some progress have also been achieved in the orthogonal and symplectic cases: Brezin-Hikami 
[TT] H8] computed the bulk scaling limit asymptotics of the moments (det fc (a — H)) and also 
provided some asymptotic expressions for averages 1)1.1. 2|) with small number of factors, and the 
asymptotics of the negative moments (l/det m (/3 — H)) has been obtained by Fyodorov-Keating 
[33] and Forrester-Keating |2"5] . 

However, the problem of computing the bulk scaling limit asymptotics of general averages 1)1.1.2)1 . 
despite considerable interest of physicists, see e.g. Andreev-Simon Gronqvist, Guhr and Kohler 
[37] . Fyodorov (3J, Szabo [51], Splttorff-Verbaarschot [52] . Zirnbauer [SHUSH], remained open. 

The main goal of this paper is to provide explicit (determinantal or pfafnan) expressions of 
averages (jl.1.2)) in terms of those that involve only one or two determinants, and to evaluate the 
asymptotics of 1)1.1. 2|) in the bulk scaling limit regime of the Gaussian ensembles in the middle of the 
spectrum. (In the case of real symmetric matrices we consider only the case when the matrices are 
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of even size.) It is worth noting that the standard correlation functions of the matrix ensembles 
can be easily extracted from averages p.l.2|) . Thus, we obtain the classical determinantal and 
pfafnan formulas for the correlation functions (see e.g. Refs. 01], [57]) as a corollary. 

1.2. The results: algebraic part. In order to state the results we need to introduce some 
notation. Take a positive measure [i on R with finite moments and infinite support. Then 

(1.2.1) C x ] '- = W\I II \xi-Xjfv{d Xl )®---®ii{dx N )^U 

for any N > 1 and (3 > 0. Take the probability measure on R^ given by 

(1.2.2) p^(dxi, . . . , dx N ) = ^—j- ] | \xi — Xjf /J,(dxi) <g) • • • <g> n{dxx). 
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N l<i<j<N 



and for any symmetric function g : R^ — > C set 

{g) A W ■= / g(x 1 ,...,x N )pP(dx 1 ,...,dx N ) 



provided that the integral converges. Also, for ( £ C set -D(C) := E[j>i(C ~~ x i)- The number of 
factors in such products will always be clear from the context. 

It is well known (see e.g. |44j ) that the radial parts (that is, projections onto different eigenvalues 
{xi}) of the probability measures 

P(dH) = const • exp(-Q(H)) dH, 

where Q(x) is an even degree polynomial with positive highest coefficient, and H belongs to the 
linear space of 2N x 2N real symmetric, N x N Hermitian, or N x N quaternion real Hermitian 
matrices 1 , are exactly the measures p^, P^\ aim Pn W1 th our measure \i being equal 



fj>(dx) 



exp(—Q(x))dx, (3 = 1,2, 
exp(-2Q(x))dx, (3 = 4. 



In the cases (3 = 1 and 2, the characteristic polynomial det(£ — H) of the random matrix H is 
exactly our product-function D{Q), while in the case (3 = 4 we have det(£ — H) = L> 2 (£). For any 
finite sets A = {ai , . . . , a p }, B = {b±, . . . , b q } denote 



H(A;B)= 11 (a, -b,), V(A) = ]] C°* ~ V ( B ) = II & ~ b ^ 

i=lj=l !<*<i<P l<*<i<9 

Theorem 1.2.1. (orthogonal and symplectic cases) (i) For any integers N > 1 and S > 1 — N, 
and finite sets of mutually distinct complex numbers 

(1.2.3) a = {ai,...,a k }, (3 = {(3%, . . . ,(3 m }, k-m = 2S, 

such that (3 n R = 0, one has 



w£\a,(3\a,(3) 



^Matrix elements of a quaternion real matrix are 2x2 matrices of the form 



where z and w are complex 



z w 
~w z 

numbers. An N x N quaternion real Hermitian matrix has 2A^ eigenvalues which come in pairs of coinciding real 
numbers. 
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where wff is a skew- symmetric (k + m) x (k + m) matrix with rows and columns parameterized 
by elements of a and (3, and with matrix elements given by 



a 



(i) 



W^'(ai,aj) = m (cti - ctj) (D(ai)D(aj)) m 



(1.2.5) W$\ai,Pj 



'2N+2S 

1 /Dlou 



at - (3j \D(f3j)/ A (i) 

lA 2N+2S 

,(1) 



(1) 



(ii) For any integers N > 1 and S > 1 — N , and finite sets of mutually distinct complex numbers 
U.2.31 such that (3 n R = 0, one has 

ri9fi x / nii^ 2 K> \ ggls n^g) pf 

where wffi is a skew- symmetric matrix with rows and columns parameterized by elements of a and 
f3, and with matrix elements given by 

(4) 



W^( a ,f3\a,P) 



\— /a (4) 



(1.2.7) ' /D2(ad 



Oi-ft \D 2 (ft)/ A M 

"iV+S 

(4) 



c 



(4) ^, \D 2 (/3i)D 2 ((3 '•) / A (4) 

Note that the formulas for the orthogonal and symplectic ensembles above are essentially iden- 
tical, except for the number of distinct eigenvalues, which in orthogonal averages is equal to twice 
that number in the symplectic averages, and for the form of the characteristic polynomial, which 
is equal to in the orthogonal case and to -D 2 (C) in the symplectic case. 

The formulas in the theorem above require that the total number k + m of factors in average 
2jl is even. These formulas can be easily extended to the situation when the total number of 
factors is odd. In order to do that, one needs to take the even case and send one of the parameters 
oti,Pj to do. After the limit transition, some of the two-point averages in the matrix elements of 

wffl will turn into one-point averages of the form (det(£ — H)), (1/ det(C — H)}. 

Interestingly enough, the analogous result in the unitary case looks more complicated! (Al- 
though, its proof is, actually, simpler.) A large part of the next result is contained in the previous 
work of Brezin-Hikami [15) I16j. Strahov-Fyodorov [SSj, Baik, Deift and Strahov j^j. 

Theorem 1.2.2. (unitary case) For any integers N > 1 and S > 1 — N, and finite sets of mutually 
distinct complex numbers 

oT = {cq,. . . ,a^J, a + = {a^, . . . ,aj^}, 

/r = {/3r,...,/?- 2 }, /?+ = {/?+,...,/?+}, 

with \a~\ - \a + \ = \/3~\ - = S, such that a + nK = 0, (3 + n R = 0, one has 

\njL^«)n?i^(/3;)/ A g) ' 

(2) 



(1.2. 



cy+s n(»';» + )n(r;/3 + ) 

C m V(a-)V(a+)V((3-)V(P+) aet 



W%\a-,f3 + \p-,a 
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is a matrix with rows parameterized by elements of a and (3 + , columns 
parameterized by elements of (3~ and a + , and with matrix elements 

,(2) 



w 



(2) 



N 



(1.2.9) 



W 



(2) 



A? 




N+S-l 



y N+S 



N+S+l 



Once again, the formula above holds for the even total number of determinants, but the odd 
case is easily obtained by sending one of the parameters af , /?■ to infinity. 

It is worth noting that Theorem 11.2.21 provides many different expressions for the same averages, 
depending on how we split the factors in the numerator and denominator into groups. The resulting 
identities are often not easy to prove independently. 

Define the nth correlation measure of pffi by 



P nN (dxi,...,dx n ) 



n 



I 



Px\dxi, . . . ,dxN)- 



For a function /(£) of a complex variable £, which is continuous in both half-planes Im£ > and 
Im£ < up to the real axis, we will denote by [/(C)]c=s> x £ the difference of the limit values 
of /(C) as C —> x from bottom and from top divided by 2-7ri: 

[f(C)] x = ^-(f(x-iO)-f(x + iO)). 
The next statement is an easy corollary of Theorems 11.2.21 and 11.2.11 

Corollary 1.2.3. Take any n > 1 and assume that near n points x%, . . . ,x n G M the measure fi 
is absolutely continuous with respect to the Lebesgue measure, and its density their is uniformly 
Holder continuous. Then the nth correlation measure of A^ (or A^y for (3 = 1) has a continuous 
density near {x\, . . . , x n ) which is given by 
• For (3 = 2 



Pnm Xl > ■■■,x n ) = det[K$(xi,Xj)] 



]JH,j=l 



where for x ^ y the kernel is given by 

K {2 \x,y) 



D(x) \ 
D(0/ 



A 



C2) 



and for x = y the kernel is defined by continuity. 
• For (3 = 1 

/ ? n,2Af( Xl ' • • • i Xn ) = ^ [-^2iV ' X j)]i,j=l 

where the skew-symmetric 2x2 matrix kernel for x ^ y is given by 



a 



m 



a 



a) 



(x-y) (D(x)D(y)) (1) 

i— in - - 



l 

x-y 



D(x) 



A 



(i) 



<=y 



l 

x-y 



c: 



(1) 



C. 



(i) 



{x - y) 



D(x)\ 

' 2N 



A 



(i) 

2JV+2. 



C=x,v=y- 
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and for x = y the kernel is defined by continuity. 
• For = A 

p%{ Xl ,...,x n ) = Vl [K^\xi,xj)]2 J=1 
where the skew-symmetric 2x2 matrix kernel for x ^ y is given by 



K%\x,y) 



°JV-1 

9r ,(4) ai \ — \m i £v 

2G N "jv+i 



{x-y)(DHx)D^y)), i) 



1 

" 2(x-y) 



D*(x)\ 



C 



(4) 
JV+1 



2C' 



(4) 



2(^-y) 



D 2 (z) 

° 2 (0/ A (i> 



C=x,v=y- 



and for x = y the kernel is defined by continuity. 



If instead of asking for the measure // to have nice density near x\, . . . , x n , we require that fi is 
purely atomic near these points, then the formulas for the correlation functions above will continue 
to hold if we understand the symbol [/(C)]as as the residue of the meromorphic function /(C) at 
the point x. 

Using standard techniques of RMT, one can evaluate the two-point (k + m = 2) averages Ql. 1.2)1 
via the (skew)-orthogonal polynomials associated with the problem, see Sections 12.71 13.51 and 13.61 
below. Then the formulas above yield a new proof of the well-known evaluation of the correlation 
functions of pffi with = 1,2,4, in terms of (skew)-orthogonal polynomials (see e.g. |56|). 

Let us emphasize that our approach provides a proof of the determinantal and pfaffian formulas 
for the correlation functions which does not use the (skew)-orthogonal polynomials. 



1.3. The results: analytic part. Theorems 11.2.11 and 11.2/2*1 are very convenient for computing 
the limits of averages p. 1.2)1 as the size of matrices goes to infinity. In what follows we will use 
the notation GOE^, GUEn, GSEn to denote the measures pffi with 



= 1, p(dx) 



2 dx, = 2, fj,(dx) 



e x2 dx, 



(3 = 4:, n{dx) = e x dx, 



respectively. The normalization is chosen in such a way that the needed scaling will be the same 
in all three cases. Note that in the case of Gaussian weights, the constants cffi are explicitly 
computed as certain products of T-functions, see e.g. Ref. [H]. The values relevant for our three 
cases will be given in Section ["3 



Theorem 1.3.1. (i) With the notation and assumptions of Theorem M. 2. lV i). we have 

rfc 



lim 



,(1) 

y 2N 



(2N)' 



2N+2S 



^/ nr=i^/v2AQ \ 

GOE(2N) 

n(a;/3) 



V(a)V(f3) 



Pf 



?GOE 



(a,P\a,P) 



where §GOE( a , P\ a , P) * s a skew- symmetric (k + m) x {k + m) matrix with rows and columns 
parameterized by elements of a and 0, and with matrix elements given by 



§>GOe(o. p , a q ) 
§>GOE{ot p , P q ) 

§>GOe(P P , Pq) ■ 



1 d sin(a p — a q ) 

7T dcti Ct n — Ct n 



a p - 

expi(f3 q -a p ) 
P q -a v ' 
expi(a p -(3q) 
a P -/3 q ' 



ImPq > 0, 
Im Pq < 0, 



' J+°° cxp(i(/3 p -/3q)t) ^ 

2~i < - j +oc cx P('(A?-ft>)*) di 



Im Pp > 0, lmP q < 0, 
Im0 p < 0, Im 0g > 0, 
in all other cases. 
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(ii) With the notation and assumptions of Theorem \1.2.1lf ii). we have 



lim 



(2NY 



N+S 



nT=iD 2 (P t /V2N)/ GSEiN) 



n(a;/9) 



Pf 



&gse(ol, P\a, P) 



V(a)V(P) 

where §gse(oi, P\a, P) is a skew- symmetric (k + m) x (k + m) matrix with rows and columns 
parameterized by elements of a and P, and with matrix elements given by 

1 f 1 sin(a p — a q )t 



§GSE{oi pi Pq) 



7T 







t 



dt, 



expi(P q -a p ) 
P q -a v 
expi(a p -l3 q ) 



$GSE(Pp,P q ) = 2m* 



Im P q > 0, 
Im P q < 0, 

Im P P > 0, Im p q < 0, 
.... W i™ r- 

q Pq~Pp 



O-p-Pq 

d expi(P P ~P q ) 

dP p P v -f3q ' 

di3 q p q -p p ' impp ^ u, impg ^ u, 
v 0, m aZZ other cases. 



Similarly, Theorem 11,2.21 yields 
Theorem 1.3.2. With the notation and assumptions of Theorem \1.2.2\ we have 



lim 

N-*oo M l 

V 



^ N {2N)~- 



N+S 



, UT=i D(a+/V2N) UT=i D(PpVW) I 



GUE(N) 



(-1) 



(|/3-| + |a-|) 2 +(|/3-|-|g-|) fJ( Q - . a +) ]"[(£- . p 



■det [§GUE(a-,P + \p-,a + )} 



V(a-)V(a+)V(0-)V(JP 

HereE>GUE(&~ , P + \P~ ,ct + ) is a matrix with rows parameterized by elements of or and P + , columns 
parameterized by elements of P~ and a + , and with matrix elements 

1 sin(a" - p~) 



^gue{u p ,P q ) 

^GUE{otp ,Oig) 
^GUE{Pp , Pq) 



K Ctr, 



Pq 



^GUEKPp 



expi(aj-a p ) 

OLq-OL p ' 

expi(a p ~ctq) 

otp— at ' 

cxpi(/3+-/3~) 

-PI ' 
expi(P q -pj) 

Pq-Pi ' 



expi(pj —aj ) 

Pi -at ' 
cxp i(aj ~pj) 

^ -Pi 



Ima| > 0, 
Ima| < 0, 



lmP+ > 0, 
Im/3+ < 0, 



2iti < 



0, 



Im/3+ > 0, Ima+ < 0, 

Im/3+ < 0, Ima+ > 0, 
in all other cases. 



Observe that if in any of the three cases (GOE, GUE, GSE) the number of determinants in 
the numerator is the same as the number of factors in the denominator, then no prefactors of 
the form /C^ or of the form 2N raised to a power are needed, the limit exists without any 
normalization. We believe that this is a deep fact, however, at the moment we do not have a 
conceptual explanation for it. 

Another intriguing observation is that if we take the same number of determinants in the numer- 
ator and denominator, and the signs of the imaginary parts of all the arguments in the denominator 
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are the same, the limit turns out to be the same in all three cases! More exactly, 
/ det(ai/\/2iV — H) ■ ■ ■ det{a k /V2N - H) \ 



lim 



\ det(/V\/2iV - H) ■ ■ ■ det(/3 k /V2N - H) / 



GOE(2N), GUE(N), GSE(N) 
k 

exp ±i ~^2(oij — (3j 
3=1 

where the positive sign inside the exponential corresponds to the case Im(3\, . . . , Im (3 k < 0, and 
the negative sign corresponds to the case Im (3\, . . . , Im (3 k > 0. Again, it would be nice to have a 
conceptual explanation for such an unexpected coincidence. In particular, one may ask if such a 
limit exists for arbitrary (3 > 0, and if so then whether it depends on (3. 

It is worth noting that if one formally uses the asymptotics of Theorems II .3. .3.21 in Corollary 
11.2.31 to compute the asymptotics of the correlation functions then one easily recovers the well- 
known results for the bulk scaling limit of the correlation kernels, see e.g. Refs. [22j, 03]. For 
details of this computation see Remark l5.1.8l Section [5. II below. 

An important feature of the bulk scaling limit in random matrix models is that the limiting 
values of the correlation kernels turn out to be universal: they depend on (3 but do not depend on 
the potential Q(x). This universality property has been verified in full generality for (3 = 2 case in 
Refs. [23J, [H]; and in the middle of the spectrum for (3 = 1, 4 in Ref. We expect that the bulk 
scaling limits of averages ljl.1.21 which we computed in Theorems 11.3.11 11.3.21 in the middle of the 
spectrum for Gaussian potentials, are also universal in a similar sense. For (3 = 2 this was proved 
in Ref. [53], and for (3 = 1,4 the question remains open. 

1.4. The method: discrete approximation. In recent years we have seen a lot of progress 
in understanding discrete probabilistic models of random matrix type which come from various 
domains of mathematics. One important observation that becomes clear from the point of view 
of those discrete models is that the measures pffi with (3 = 1,2,4, may be viewed as degenerate 
cases of more general objects, the so-called determinantal point processes in the unitary case and 
pfaffian point processes in the orthogonal and symplectic case. 

Inside these more general classes of point processes one also finds discrete analogs of the measures 
pffl- In the unitary case, the discrete ensemble is different from the continuous one only by the 
fact that the support of the measure fi is discrete. In the orthogonal and symplectic cases the 
difference is more substantial. 

Let X be a finite subset of M with an even number of elements. The discrete analog of the 
measure p^jv ^ s a measure on 2iV-point subsets X of X which have the following property: for any 
x £ X the number of points in X that are smaller than x has the same parity as the number of 
point in X that are smaller than x. The weight of one such subset is given by the familiar formula, 
cf. (TT2~2|l 

2N 

(1.4.1) Prob(xi, ... ,x 2 n) = const }~[ \xi - Xj\ ■ JJ /(xj), 

l<i<j<2N i=l 

where / is some positive weight function. We call this measure a discrete orthogonal ensemble. 

A discrete symplectic ensemble is defined by the same formula (|1.4.1|) . but the admissible subsets 
X are different: with any point x such subset must contain the immediate predecessor of x in X. 

It is not hard to see that if we take X to be a lattice then in the limit when the step of this 
lattice goes to zero, the discrete orthogonal/symplectic ensembles turn into measures similar to 
j»2jv an d P$ • (Indeed, if the points of X are split into pairs, and in each pair the points are 
innnitesimaily close, then the Vandermonde determinant of size |X| in (jl.4.1|) is asymptotically 
given by the fourth power of the Vandermonde determinant of size |AT|/2 in scaled locations of the 
pairs.) 
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The discrete orthogonal and symplectic ensembles are related by the particle-hole involution: if 
X is from a discrete orthogonal ensemble than X = X \ X is from a discrete symplectic ensemble 
(with a different weight function and different number of particles though) and vice versa. The 

(2) 

particle-hole involution of a discrete unitary ensemble (which is a measure of the form p v N ' with 
discretely supported /i) is again a discrete unitary ensemble. In the (3 = 2 case the particle-hole 
involution has been used before, see Refs. J2] (Section 5), [Hj, [22], 

The motivation for introducing these discrete objects comes from representation theory; in cer- 
tain models these are just different ways of parameterizing partitions. If A = (Ai,A2,...) is a 
partition then {2Aj — i} is an admissible point configuration for a discrete orthogonal ensemble, 
while {Aj — 2i, Aj — 2i + 1} is an admissible point configuration for a discrete symplectic ensemble. 
The convenience of coordinates {Aj — (3i/2} in the models involving Vandermonde determinant 
raised to the power (3 can be observed, e.g., in Refs. [33], QUI- The fact that discrete orthogonal 
and symplectic ensembles are dual to each other comes from two different ways of parameterizing 
partitions: using rows or columns of the corresponding Young diagram. We hope to give more 
details on this connection in a future publication. 

In the paper we actually prove analogs of Theorems 11.2.11 and 11.2.21 for discrete orthogonal, 
unitary, and symplectic ensembles and then obtain Theorems 11.2.11 and 11.2.21 by limit transition. 
Note that two parts of Theorem 11.2.11 in the discrete situation are simply equivalent — they are 
obtained one from the other by the particle-hole involution. This explains the similarity of the 
formulas. 

The proofs in the discrete case are obtained as follows. It turns out that if we apply the particle- 
hole involution on a suitable part of the phase space X then the orthogonal/symplectic and unitary 
ensembles turn into what we call L~ensembles. This is another subclass of pfaffian/determinantal 
point processes; the processes in this subclass are defined using a matrix of size \X\ x \X\ which is 
often denoted by L. The discrete analogs of Theorems 11.2.11 and 1 1 . 2 . 21 inter preted in the language 
of this matrix L provide some relations between matrix elements of L and its resolvent. These 
relations are then proved using standard linear algebraic facts. 

We would like to emphasize that after the objects of interest are represented in terms of L- 
ensembles, the proofs become very simple. The nontrivial part of our approach is in constructing 
the discrete models which make the statements easy to prove. 

1.5. The asymptotics. Despite the fact that Theorems 11.3.11 and 11.3.21 look like rather simple- 
minded corollaries of Theorems 11.2.11 and 11.2.21 rigorous proofs of the asymptotic formulas do 
require some efforts. 

The first step of our computation consists of expressing the two-point averages - matrix elements 
of wffl in Theorems 11.2.11 and 11.2.21 — in terms of associated (skew)-orthogonal polynomials. 
As is well known, in all three Gaussian ensembles the corresponding polynomials are expressible 
through classical Hermite polynomials, see e.g. |44[ I27 | IT] I47j. Also, certain summation formulas 
of Christoffel-Darboux type for all three case (3 = 1,2, 4 are known, see 

EE]- Using these 

formulas and applying some algebraic transformations (which are least simple in the symplectic 
case), we write the two-point averages as finite expressions involving Hermite polynomials, their 
derivatives, and integrals. 

Further computations are easy to explain — we replace the Hermite polynomials of large degree 
by their asymptotic oscillatory behavior near the origin. However, sometimes we need to do that in 
double or triple integrals, and the resulting limit integrals are often only conditionally convergent. 
Such actions require accurate estimates of the error terms, and they constitute a large part of 
Section [5] In these estimates it is essential for us to have uniform asymptotics of the Hermite 
polynomials on the real line. The needed results follow from much more general results of [231 121] 
2 ; we adopt their formulas to the Hermite case. 



'The classical estimates proved in |55| turn to be to be insufficient for our purposes. 
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1.6. Final remarks. To conclude the introduction we would like to suggest the following debatable 
point of view which we developed while working on the subject: 

The averages of products and ratios of characteristic polynomials are more fundamental charac- 
teristics of random matrix models than the correlation functions. 

This is especially visible in the (3 = 1,4 cases. Here are some arguments to defend this thesis: 

• The pfaffian/determinantal formulas for these averages are simpler than those for the correlation 
functions. The formulas express many-point averages in terms of one- and two-point ones, and 
they do not require the introduction of (skew)-orthogonal polynomials or Christoffel-Darboux type 
kernels. 

• The formulas for the correlation functions are easily recovered from those for averages (1). 

• Even though there is no clear probabilistic sense in taking the bulk scaling limit of such averages, it 
does exist after proper normalization, and it appears to be fairly universal. The limiting correlation 
kernels are easily obtained from these limits. 

We hope to provide further arguments related to discrete probabilistic models arising in repre- 
sentation theory in future publications. 

1.7. Acknowledgements. We are very grateful to Grigori Olshanski for numerous discussions of 
the (3 = 2 case. This research was partially conducted during the period one of the authors (A.B.) 
served as a Clay Mathematics Institute Research Fellow. 

2. Determinantal point ensembles 

2.1. Point configurations. Let X be a finite set. We will denote by X,Y, ... subsets of the set 
X and we will call them "point configurations". Let X have a fixed splitting into the union of two 
disjoint subsets ("positive" and "negative"), 

X = X_ U X+ 

Then any point configuration X is a unit of two disjoint sets as well, 

X = x_ux + 
x_=xnx_, x + = xnx+ 

and we will say that the point configuration X consists of positive particles (elements of X + ) and 
negative particles (elements of X—). If a given configuration includes an equal number of positive 
and negative particles we will say that such a configuration is a balanced point configuration. Thus 

\X+\ = \X_\ 

if X = (X-\X + ) is a balanced configuration. Here \A\ denotes the number of elements in the set 
A. 

If X is a point configuration of particles in X its complement X \ X in the set X will be a 
point configuration as well. It is natural to refer to the point configuration X \ X as the point 
configuration of holes. The set (X \ X) n X_ will be called the set of positive holes, and the set 
(X \ X) n X-|_ will be called the set of negative holes. 

To any configuration X of particles there will correspond (in many ways) a balanced point 
configuration Z, which consists of both particles and holes. Namely let X = (X-\X + ) be the 
splitting of the point configuration X into positive and negative particles, and assume that |-X"_| — 
\X+\ = S > 0. Then the corresponding balanced configuration can be constructed in accordance 
with the formulas: 

z = z_ u z+, z + = z n (x+ u x ), z_ = z n (x_ \ x ) 

Z— = X— \ Xo 

Z + = X + U (Xo \ X) 

where Xo is an arbitrary subset of the set X_ consisting of S points. Fig. 1 shows the decomposition 
of the set X, where the set Xq is chosen to be the right-hand subset of the set X_. Fig. 2 
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represents the unbalanced particle-particle configuration X, and Fig. 3 explains the construction 
of the balanced particle-hole configuration Z. 

It is not hard to see that the point configuration Z constructed in such a way is a balanced 
point configuration, i.e. |Z+| = \ZJ\. Indeed if |A + | = d then |-X_| = S + d. Assume that the 
chosen subset Xq includes d\ negative particles of the configuration X. It means that Xq includes 
S — d\ positive holes. The set Z + consists of these positive holes and positive particles of the point 
configuration X, thus \Z + \ = S — d\ + d. Since Z_ consists of the negative particles lying outside 
the set Xq , \ZJ\ = S + d — d\. Therefore | Z + \ = \Z_\, and the configuration Z defined above is a 
balanced configuration. 

We will say that the balanced configuration Z is obtained from the configuration X by the 
particle-hole involution on the set Xq. 

2.2. Definition of L-ensembles. First we recall the notion of determinantal L-ensembles follow- 
ing Daley and Vere- Jones Borodin and Olshanski 112j . A random point process on the 
space X is defined by an introduction of a probability Probi(X) for each subset X of X, so that 

Prob L (X) = 1, Pvob L (X) > 

xc x 

Let L be a X x X matrix whose rows and columns are parameterized by the points of X. We use 
the notation A(a\f3) for a matrix A and subsets a = (a%, . . . ,a a ), (3 = (/3i, . . . , (3b) of its rows 
and columns, to denote the submatrix ||^4 Q .^.|| of A. Then to any subset X (or for any point 
configuration) there will correspond the diagonal minor det L(X\X) of the matrix L. Assume 
further that any such diagonal minor of L is nonnegative. 

Definition 2.2.1. A random point process living on the space of point configurations in X is called 
an L-ensemble, if 

It is clear that Yl P roD L(A) = 1 as the normalization constant det(l + L) has the following 

xc x 

well-known decomposition: 

det (1 + L) = ^2 det L{X\X). 
xc x 

According to the decomposition of the set X into the union of positive and negative subspaces, 
X = X- U X+ we write the matrix L in the block form: 

L+- L ++ 

Let h be a real function defined on X. We are interested in L-ensembles for which the matrix L is 
given by the formula 

A 



(2.2.1) L 



(2.2.2) L 



-A T 



A( X ,y) - Hx)h{y) 



x-y 

In that case we use the formula for Cauchy determinants to rewrite the expression for the probability 
of a random point configuration as 

f i v 2 (x + )vHx^ _ 

(2.2.3) Prob L (A) = <^ det(l + L) \\ 2 {X^X + ) 1 h 1+1 ' 

[ 0, otherwise, 

k I 

where n(^' B) = f] FJ (a^ — bj) for any two sets A = (ai, . . . , a*.), B = (b\, . . . , 6;), V{X) is the 
i=ij=i 

Vandermonde determinant associated with the set X, 

V{X)= JJ (xi-xj), X = (xt, . . . ,x N ) 

l<i<j<N 
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N 

and h(X) = f\ h(xj). 

2.3. L-ensemble. We have seen in the section above that an L-ensemble is completely determined 
by a splitting of the discrete set X into positive and negative parts, and by a weight h on X. With 

a different splitting of the set X, X = (^X-\X-^j and different weight h we will construct another 

ensemble, which we will call L-ensemble. 

Let X- = X- \ Xo and X + = Xq U X +1 \Xq\ = S. As we have already seen in Section |2~T1 to the 
unbalanced configurations X of particles 

X = (X-\X+), - \X+\ = S 

with respect to the splitting X = (SE_|3f+) there corresponds the balanced hole-particle configura- 
tion Z, 

Z = (Z\Z+), \Z\ = \ZA 



with respect to the new splitting X = ^3t_|X+J . The point configurations X and Z are related by 

the hole-particle involution on the set Xo as it is described in Section 12.11 

Consider an L-ensemble on X with the matrix L having the following block structure: 



(2.3.1) 



L 





-A T 



A 




A(x,y) 



h(x)h(y) 



y 



Here the decomposition of the matrix L into the blocks corresponds to the new splitting X = 
(^X-\X + ^j of the discrete set X. As soon as the particle-hole configurations Z are balanced with 
respect to this splitting we find 



(2.3.2) 



Prob,(Z) 



1 



h\Z) 



V 2 (Z_)V 2 (Z^ 



det(l + L) ^ ' \Y{Z ;Z^ 

The balanced configurations Z of particles and holes were constructed from the unbalanced config- 
urations X of particles, in which the number of negative particles is larger by 5 then the number 
of positive particles. Clearly, it is possible to rewrite Prob^(Z) in terms of the configurations X. 
With a suitable choice of the weight h the expression for Prob^(Z) in terms of configurations X 
takes the same form as the right-hand side of equation (|2.2.3jl (up to a normalization constant). 
Namely, introduce the new weight h in terms of the old weight h according to the formula: 

y), zeX-\X , 



(2.3.3) 



h(z) 



Hz) n (*- 

ye Xo 

1 



h(z) 



z e X , 



z G X + . 



I Uyex ( z -y) 

Proposition 2.3.1. This choice of the weight h (equation \2.S. 1$ ) gives 

V 2 {X_)V 2 {X^ 



(2.3.4) 



VmhAZ) 



1 



1 



h\X) 



LK "' det(l + L) V 2 (X )h 2 (X ) " X\ 2 (X_ ] X + , 
Proof. In order to see that equation (|2.3,4j) is valid we rewrite the expression 



h\X) 



V 2 (X_)V 2 (X^ 
U 2 (^;X + ) 



in terms of the balanced configuration Z (the relation between the unbalanced configuration X 



and the balanced configuration is shown in Tab. 1. We will also use the notation Z+ 11 introduced 
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there). In particular we find 

. v 2 (x^) = v 2 (x \ z{) ■ v 2 (z_) ■ n 2 (Z-;*o \ z{) 
. V2{Z+) 



v\z\)-v?{zl-z\ l ) 
. n 2 (*-; x+) = Yfiz^zl 1 ) ■ n 2 (*o \ Z{; z?) 

Thus 

V 2 {X^)V 2 (X + ) _ V 2 (Z_)V 2 (Z + ) v 2 (x \z{).u 2 (z^x ) 



U 2 (X-;X + ) U 2 (Z-;Z + ) V 2 (Z[) ■ U 2 (Zi; Z%) • Tf^o \ Z{; Z?) 

We rewrite V 2 (X \ Zl.) 



V\X \Zl) = —— ^ 



V 2 (Zl).il 2 (Zl;X \Z[) 
and note that the resulting denominator equals 

V\ZD ■ Z») ■ n 2 (X \ Zl; Z{) • jfC*) \ Z r +] Z?) ■ V 2 {Z{) 

= V 2 (Zi) • n 2 (^; 3Eo) • n'^o \ Zl; Z{) • V 2 {Z{) 

= n {x-yf- n (x- y ) 2 - n {x-y?- n ( x -y) 2 = n ^-y) 2 - 

x, ye z' + xe z' + xe z[ T x, ye z\ xe z+ 

y<x y£ x \Z[ ye Xa V> x y& v ^ x 

Introducing h(z) by equation ()2.3.3jl we see that equation (|2.3.4|) is equivalent to equation (|2.3,2[) . 

□ 

2.4. Correlation functions for L-ensembles and averages of characteristic polynomi- 
als. By correlation functions for L-ensembles g(X) we mean the probabilities that the random 
configurations include fixed sets X. Thus 

e (x) = £ Prob L (y). 

YDX 

It is known that the correlation functions g(X) are expressed as symmetric minors of the correlation 
kernel K = L(l + L)" 1 , i.e. 

(2.4.1) g{X) = det K{X\X). 

The proof of this fact can be found in the book by Daley and Vere- Jones |21| (Exercise 5.4.7), 
and in Borodin and Olshanski (Proposition 2.1), Borodin, Okounkov and Olshanski ^Uj (Ap- 
pendix). Equation (|2.4.1|) tell us that random processes associated with L-ensembles are discrete 
determinantal point processes (see Borodin and Olshanski [111112] for definitions. A comprehensive 
survey on determinantal point processes is given in Soshnikov |51j). 

In this Section we express the minors of the matrix K in terms of averages of "characteristic 
polynomials" E(a, X) associated with L-ensembles. These objects were first introduced in Borodin 
and Olshanski |13j and are constructed as follows. With the above notations, for any fixed set a 
and point configurations X and Z we set 

E(a X) - E ^ Z) 
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In order to compute the minors of the matrix K we need the following Lemma. 
Lemma 2.4.1. Let A, B be two balanced configurations on the set X, 



\A. 



Then 



where 



det [L(A-,B+\B-,A+)] = (- 



w a b = ab + b + 



| = a, \B + \ = 

w " b h{A)h{B) 
a(a 



1) 



\B-\ = b 

V(A-.)V(A + )V(B-)V(B + ; 
U(A-;A+) U(B-;B + ) 

-1) 



+ 



Proof. The formula above follows from the explicit structure of the matrix L (equation 1)2.2.2)1 ). 
from formulas for the Cauchy determinants (equation 1)6.1.1)) ') and the determinants of the block 
matrices (equation 1)6.2. 2JI ). □ 

Introduce nonintersecting sets a^, (3^ of complex numbers with nonequal elements, 



o 



(2.4.3) ct 

(2.4.4) a" = (of,... ,a mi ) 
Assume that n X = 0, n X = 0, and 

\a + \ - \a~\ = \(3 + \ - \p-\ = S. 

In what follows we extend the definitions of the matrices K, L to the sets , (3^ in the following 
way. We add to L rows parameterized by a~ U (3 + and columns parameterized by a + U (3~ , and 
then define new matrix elements of L according to 1)2.2.1)1 - 1)2.2.2)1 . where we assume that a~ and 
f3~ are added to X_, and a + , j3 + are added to X + . Then we set extended K to be related to the 
extended L by K = L(l + L)~ x . 

Proposition 2.4.2. The minor det [K(a~, f3 + \(3~, a + )] can be given as a normalized average of 
a ratio of the functions E(.,.) introduced above with respect to the L- ensemble. Namely, 



(2.4.5) 



det [K{a-,/3 + \(3-,a + )] = (-) w ^ 

h(a 



det(1 + i) h\x o) vHx 



det (1 + L) 
V(a~)V(a+) 



n(« ;« 

n(«-;Xo)n(/3-;X ) 

n(a+;X )n(/3 + ;^o) 
a" I]. 

Proof. Equation 1)6.5.2)1 gives 
(2.4.6) det [K(a~, f3 + \(3~, a + )] = ^ 



h(P) 



V{(3-)V{(3 + 

E(a+,Z)E(0+,Z) 
E{a~,Z)E((3-,Z)/ L 



where w a ,p = ~[(\a \ + |/3 |) 2 + \/3 



det L(a~,l3 + ,X\p-,a+,X) 



xc x 



We define the following sets 
(2.4.7) 



^4 



det(l + L) 
a + Ul + , A_=a~UX_, 



(2.4.8) B + = {3 + UX + , 5_=^UI_. 

Only those configurations X will contribute to the sum in equation 1)2.4.6)) for which the sets A_ , 
A + have equal number of elements and the sets B_ , B + have equal number of elements. Otherwise 



det [L{a-,p + ,X\(3-,a + ,X)] = 0. 
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It means that the sum in equation q2.4.6JI runs over unbalanced configurations X, 

IX J - \X . I = s. 



Now we apply Lemma 12.4.11 to compute det [L{a~, (3 + , X\j3~, a + , X)] with A±, B± given by 
equations (|2.4.7(l and (|2.4.8j) respectively. We find 

-i 



det [K(a-,f3 + \f3-,a> 



(_)»-,/> [det(l + L)] 
V(a-)V(a+)' 



U{ a ~> a+ ) 
E(a + ,X)E(P+,X 



h(/3) 



V(f3-)V{{3 + ) 

n(/?-;/3 + ) 



E(a-,X)E(p-,X) 



/* 2 (X) 



V(JT-)V(A' + ) 



n 2 (x 



Since the sum in this equation runs over unbalanced configurations, we cannot interpret the sum 
as an average over the L-ensemble. However, it is possible to rewrite this sum in terms of the 
balanced configurations Z of particles and holes as it was explained in Section \'2.'dl Then equation 
Q2.3.4JI tell us that the sum in equation above is just an average over L-ensemble. Replacing the 
sum over unbalanced configurations X by the sum over balanced configurations Z and using (J2.4.2[) 
we prove the Proposition. □ 

Remark 2.4.3. The Proposition above clearly gives the matrix elements of the matrix K in terms 
of averages of E' s . For example, 



det (1 + L) , 2 



det (1 + L) 



h\x ) 



h{(3t)h{a\ 



{Pi ~ x ){aj - s ) 



E(f3+,Z)E(a+,Z) /L 



where the set Xq in the definition of the ensemble L (see Tab.l) consists of only one point xq, 
Xo = {xq}. The point xq is an arbitrary point in X_. The weight h of this L-ensemble is 

fh{z){z-x ), zeL\{i} 



(2.4.9) 



Hz) 



i 



h{x ) ' 

Hz) 
{z - s ) 



Z = Xq 



z £ 3C-i 



as it follows from expression (|2.3.3[) . 



2.5. Discrete polynomial ensembles. Assume that a nonnegative function /(s) is given on a 
finite set X C R. We also require that / does not vanish at least at N distinct points. Then 
the monic orthogonal polynomials ttq = 1,7Tl, . . . ,vrjv_i can be introduced by orthogonalizing the 
system (l,s, . . . jS^ -1 ) in the Hilbert space L 2 (3£,f[i). Here /j, denotes the counting measure on 
the set X. The orthogonality condition for the monic orthogonal polynomials on X is 



c? 6 



'jk- 



(2.5.1) J2 *j(x)*k(x)f(x) 

xe X 

Let Conf]\r{X) denote the set of iV-point configurations (subsets) in X. Such configurations will 
be denoted by X A , X A G Confjq(?C). Consider a point process on X which lives on Conf^{X) 
and for which the probability of a configuration X is given by 

(2.5.2) Prob(X A ) = const | | f{x) • V 2 {X A ), X A eConf N (X). 
We will denote this process by Ajv(/)- 

Definition 2.5.1. The point process Ajv(/) will be called iV-point {(3 = 2) discrete polynomial 
ensemble with the weight function /. (The value of (5 refers to the power of V{X) in (|2.5.2|) .) 
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For symmetric functions g = g yx A , . . . , x^J of points of the configuration X A , the aver- 

e ( g)A N (f) W1 th respect to Prob(X A ) is 

E g{X^)f{X^).V\X^) 



, > _ X A e Conf N (X) 



E /(^ A ) • V 2 {X A ) 

X A e Conf N (X) 

Define the characteristic polynomial cZ(£) associated with the point configuration X A and the 
complex parameter £: 

(2.5.3) d(0=H&X A ). 

Then the discrete variant of the classical Heine identity is given by the following 
Proposition 2.5.2. 

n N (0 = (d(0)A N (f) 

where (-)ajv(/) denotes the average over the discrete polynomial ensemble Ajv(/)- 

Proof. The proof of this relation is an application of standard arguments of the theory of orthogonal 
polynomials (see, for example, Deift |22| . Chapter 3). □ 

Proposition 2.5.3. The m-point correlation function 

Q m (yi, . . . ,y m ) = ^2 Prob(X A ) 

X A e Conf N (X) 
X A 9 {yi,...,y m } 

for the discrete polynomial ensemble is given by the determinantal formula 

(2.5.4) g m (y 1 ,...,y m ) = det[K^ D (y i ,y j )]^. =1 , m = 1, 2, . . . 

where K^ D (x,y) stands for the normalized Christoffel-Darboux kernel, 

One possible proof of this proposition is to adopt standard arguments of the random matrix 
theory to the case of the discrete ensemble Atv(/) (see, for example, Mehta |44| . Deift jUj). In 
the following sections we present another proof. First, we will show that the discrete polynomial 
ensembles are equivalent to the L-ensembles (see Section [2~H|l . Specifically, the discrete polyno- 
mial ensembles and the L-ensembles can be transformed to each other by a suitable particle-hole 
involution. Second, we compute the averages of products and ratios of characteristic polynomials 
E(a,X) for the L-ensembles (see Section l2~7|) . From these averages we deduce the averages of 
products and ratios of characteristic polynomials d(£) for the discrete orthogonal ensemble A]y(f). 
Third, we show in Section f2.8l how formula (|2.5.4|) can be obtained from the resulting expressions. 

2.6. Equivalence of L-ensembles and discrete polynomial ensembles. The connection be- 
tween the discrete polynomial ensembles Ajy(/) and L-ensembles was demonstrated in Borodin 
and Olshanski |12j . §5 . In the terminology of Section 12.11 this relation is a consequence of the 
particle-hole involution on an TV-point subset of X. To be more specific, assume that the subset 
X+ in the definition of the L-ensemble (see Section 12.2(1 is finite, and consists of N points. Thus 
|3E+| = N. We assume that the weight h of the L-ensemble is nonnegative on X and strictly positive 
on X+. If the L-ensemble is defined by equation (J2.2.2[) then only balanced configurations X with 
\X + \ = \X^ \ have non-zero probabilities. Consider the particle-hole configuration X^ constructed 
by the particle-hole involution on the set X+. Thus 

X A = X- U (£+ \ X+) . 
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Clearly, the configuration consists of precisely N — \X+ \ + |X_ | = N points where N — \X+\ is 
the number of negative holes, and \X-\ is the number of negative particles. Moreover, it is always 
possible to introduce a weight / on X in such a way that Piob(X^) takes the same form as in the 
definition of the discrete polynomial ensemble (equation (|2.5.2jl ). Namely, introduce the weight / 
in terms of the weight h of the L-ensemble by the formula 

t h 2 {x) 



(2.6.1) fix) 



x <E X. 

X £ %, 

> h2 ( x ) n. y e3L+,vj£z( x -v) 2 ' 



Proposition 2.6.1. The discrete polynomial ensemble with the weight f and the L -ensemble 
defined by equation \2.2.2]) - \2.2. are connected by the particle-hole involution on the set X+. 

Proof. See Borodin and Olshanski ^2]) Proposition 5.2. □ 

Once the L-ensemble is given, we can construct an L-ensemble with respect to the new splitting 
of the set X, X = as in Section |2~31 (see also Tab. 1). Clearly \X+\ = N + S. This L- 

ensemble induces a discrete (N + 5)-point polynomial ensemble under the particle-hole involution 
on the set 3E+. It follows from Proposition 12.6.11 that the induced polynomial ensemble with the 
weight f(x) will be equivalent to the L-ensemble with the weight h{x) if 

h 2 (x) 



(2.6.2) f{x) 



rL P r>-y) 2 



x £ X- 



x G X-i 



Denote by A^r(/) the discrete iV-point polynomial ensemble which is equivalent to the L-ensemble, 
and by A^+sif) the discrete (N + <5)-point polynomial ensemble which is equivalent to the L- 
ensemble. 

Proposition 2.6.2. The ensembles Ajv(/) and Ajy + s(f) have the same weight, 

f{x) = f(x), x ex. 

Proof. We use formula (|2.H.H|) to express the weight /i(x) in terms of the weight h{x) in expression 
(|2.6.2j) . and after that formula (|2.6.1|) which expresses the weight f{x) in terms of the weight h(x): 
• x £ X- \ Xq 

f = h 2 {x) = h 2 (x)Y[ yeXo {x-yf 

nx) n ye x + ( x -y) 2 rw>-y) 2 

h 2 ix) 



• x £ X 



o 



II, v) 2 f{x 

f{x) _ 1 _h 2 {x)X[ yeXQ ^ x {x-yf 



h 2 {x) 



x 6 X> 



U y ex+( x -yy 
Uyex ( x -y) 2 



H x ) n y6 w > - y) 2 h2{x) n ^ w> - y)2 

= h 2 {x)\{ yex+ ^ x (x-y) 2 = /(X) 
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□ 

The correspondence between the L-ensembles and the discrete orthogonal polynomial ensembles 
is summarized in Tab. 2. 

2.7. Averages of characteristic polynomials. The goal of this Section is to compute the aver- 
ages of characteristic polynomials for the discrete polynomial ensembles. We will use the equiva- 
lence of the L-ensembles and the discrete polynomial ensembles, and deduce from Proposition [2^21 
the averages of products and ratios of characteristic polynomials. Specifically, we will express the 
right-hand side of equation (|2.4.5j) in terms of averages and normalization constants of the ensemble 
Ajv(/). This ensemble has the weight / given by equation (|2.(j.lj) . Thus A^(f) is connected with 
the L-ensemble by the particle- hole involution on the set X+ in accordance with Proposition 12 . 6. II 
We begin with the constant in equation (|2.4.5j) . Denote by Cm and Cn+s the normalization 
constants for Ajy(f) and Aat + s(/): 

C N = £ V 2 (X)-f(X), 

le Conf N (X) 

C N+ s = E y2 ( X ) • /(*)■ 

X<= Conf N+s (X) 

Proposition 2.7.1. The constant in equation \2.4-5\j is equal to the ratio of the normalization 
constants Cjv and Cn+s, i-e- 

C N+S _ V 2 (r )h 2 (r \ det (1 + L) 
~Cn~~ {X ° )h (Xo) det(l + L)- 

Proof. The expression [det (1 + L)] -1 represents the probability that a random configuration X of 
the L-ensemble is empty, i.e. 

Prob {X = 0) = [det (1 + L)Y l . 

The empty configuration (of particles) for the L-ensemble corresponds to the configuration X A = 
X+ of the ./V positive holes for the Ajv(/) ensemble. Since the L-ensemble is equivalent to the 
Ajv(/) ensemble we have 

Prob (X = 0) = Prob (x A = X+) = — V 2 {X + )f{X + ). 

Therefore 

[det (l + L)n 1 = ^-l/ 2 (X+)/(£ + ). 

We repeat the above considerations for the L-ensemble and conclude that 

det (1 + L) _ C N +s V 2 (X + )f(X+) 
det (1 + L) C N V 2 (X U X+) f {X U X + ) 
Cn+s 1 



C N F2 (Xo)n 2 (Xo;X+)/(Xo) 

Cn+s 1 



C N V 2 (X )h 2 (X )- 

Here we have used formula ()2.6.1|) to obtain the last equation. □ 

Let the sets cr* 1 , (3^ be defined as in Section l2~4l Recall that (f(£) was defined in equation (|2.5.3|) , 

Theorem 2.7.2. For any integer N > 1, take an integer S such that N —\X\ + \ < S < N—l, com- 
plex numbers aT = (aj , . . . , a'J, q+ = (af, . . . , a£), (3~ = (#f, . . . , f3~ 2 ), /?+ = (f3f , . . . , (3^ ) 
such that \a~ \ — \a + \ = \ — = S, in each set , (3^ the numbers are pairwise distinct, 
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and the sets a + ,/3 + do not intersect X. Then the average of products and ratios of characteristic 
polynomials with respect to the discrete polynomial ensemble is given by the formula 



(2.7.1) 



^■Lid(apUj 2 =id(P] 



Ajv(/) 



Cat- 



Cat 



n(« 



; a 



V(a-)V(a+)V(p-)V((3+) 
where the kernel function Wn is defined by 

Cat- 5-1 



det [W N {a-,(3 + \f3-,a + )] 



• W N (a-,f3-) 




'N-S 



and S, w a/ 3 are expressed in terms of |a ± |, \(3 ± \ 
any integer, positive or negative. 



A N _ s+1 (f) 

as in Proposition \2.4-^ Note that S above can be 



Proof. Assume first that S > 0. Then the proof is based on the expression for the minors of the 
K matrix in terms of averages of E' s over the L-ensemble (see Proposition 12. 4. 2|) . The equivalence 
of the L-ensembles and orthogonal polynomial ensembles enables us to rewrite averages of E' s as 
averages over characteristic polynomials with respect to orthogonal polynomial ensembles. We 
note that 

py 7 , = U(^z + ) _ n(Q;*o)n(«;*+) 

where Z A is a random point configuration for the ensemble Aat+s(/)- We insert the expression 
above for E(a, Z) to the formula for the minors of the matrix K in Proposition 12.4,21 According 
to Proposition I2.6.T1 and Proposition 12.6.21 the average in the formula for det [K(a~ , (3 + \f3~ ,a + )] 
can be understood as that over the discrete polynomial ensemble Aat+s(/)- Thus the formula for 
the minors of the matrix K takes the form 



det [K(a-,p+\p-,a* 



(2.7.2) 



(_)™a,/3 

n(a 



Cat+5 



Cat 

x+)n(/3 + ;X 



V(g-)V(a + ) 
Y\{a-;a+) 



n(a";X + )n(/3-;X 4 



UT=\d(a-)UT=\d(Pn \ 
. 11*=! d(ap UjLi d(P+) I AN+s{f) 



It is clear that the formula above expresses the matrix elements of the matrix K as averages of 
corresponding products and ratios of characteristic polynomials. In particular we find 



K(J3+\a+) = h(a+)h(P i 
Cjy+i / 

Cat \d(p+)d(a 




JV+1 



(/) 



K{aT\(3j) 



Cat 

x {d(ar)d((3 j 



- h(ar)h((37)l[ \ a -;X + )H \(3-;X. 



Aiv-i(/) 
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1 / d ( a i)\ 



or -a] \d(a+)/ AN{f) 



v 1 /d(pr 




gf - B~ 

^ \ 'I A N (f) 

Expressing the elements of the matrix K in the left-hand side of equation (|2.7.2|) in accordance 
with these formulas we prove the Theorem after some simplifications. 

If 5 < we can prove Proposition 12.4.21 and Theorem 12.7.21 in a similar way considering Xo in 
the construction of balanced particle-hole configurations as a subset of 3£+. Another possibility to 
see that formula (|2.7.1j) for S < holds is to perform the particle-hole involution on the whole X. 
Under this particle-hole involution formula (|2.7.1|) does note change. However, S > is transformed 
to S < 0. □ 

Remark 2.7.3. In formula (|2.7,1|) the number of characteristic polynomials is greater or smaller 
than that in the numerator by 2|5|. The case when these numbers are different by an odd integer 
is obtained from (|2.7,1|) by taking the limit aj -> oo or 0j — * oo for some % or j, and using relation 

d(e)=^(i + o(i)),^oo. 

Remark 2.7.4. Theorem l2.7.2l combines the different formulas of two-point function type obtained 
in Brezin and Hikami JH]) Strahov and Fyodorov [^Sl into one formula (equation (12.7.11) 1. 

Proposition 2.7.5. The kernels Wjy in the theorem above can be expressed in terms of the monic 
discrete orthogonal polynomials 7T/t(x) and their Cauchy type transforms /ifc(e), 

, , n 1 ftk(%) f(x) 
2.7.3 h k (e) = — 1 



2m 



In particular we find 

W N (a-,Pj 



Cn-s a i - 0j 



w / - +x n N - S (a^)h N -s-i(a i )-'K N -. S - 1 (a^)h N -s(a i ) 
WjvK ,aT) = jn-S-1 



a 



k 



yVN{Pi ,Pj ) = JN-S-1 



v , „ , hM^(3 + )h 



3 

Here we have introduced the coefficients jk, 

2m 
c k 

Proof. The computations for the case of the discrete orthogonal ensemble Ajy(/) are similar to 
that for the case of the ensemble of Hermitian matrices (see Strahov and Fyodorov |53| . Baik, Deift 
and Strahov pj). □ 
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2.8. Correlation functions of discrete polynomial ensemble. This Section shows the relation 
between m-point correlation function Q m (yi, ■ ■ ■ , y m ) for the discrete polynomial ensembles and 
averages of characteristic polynomials. The Proposition below gives Q m (yx, ■ ■ ■ , y m ) in terms of 
these averages. 



Proposition 2.8.1. For any y±, . . . ,y m £ X we have 



(2.8.1) Q m (yi, ■ ■ ■ ,Vm) = Res . . . Res 

vi=yi v m =y m 



0" 



I d{u\) . . . d{i 



dux-.-dum \ d(vi) . . . d{v m ) 



Aat(/) 



Proof. In order to see that that the expression for g m (yi, ■ ■ ■ ,y m ) is correct we differentiate the 
product of characteristic polynomials in the numerator, 



d" 



du\ . . . du Tl 



m N 



nriK-f] 

i=ii=i 



N 

£ 

jl,— ,jm=l 



An 



m N 



Thus the formula in the proposition is equivalent to the following one: 

r N 



Qm(yi, ■ ■ ■ ,y m ) = Res . . . Res 

ui=2/i v m =y rr 



^ ^ Piob(JT) ^ ] r~r m ( A \ 

Xd Conf N (X) L ji,...,j m =l llk=n v k ~~ X j k > 



which is evidently equivalent to the definition of g m . 

Remark 2.8.2. The above argument applies to any point process on X, not just Ajy(/). 



□ 



Proposition 2.8.3. The m-point correlation function of the discrete ensemble Ajv(/) is given by 
the formula 



(2.8.2) 

where 

(2.8.3) 



Q m (yi,---,ym) = det [if(yi,yj)]™,- =1 



K(x,y) = < 



Res 



Res 



d{x)\ 

±/M\ 

dr)\d(0/ 



Aiv(/) 



x^y\ 



x = y. 



Proof. We note first that if we take 

|/?+| = |/T| = 0, \a + \ = \oT\ = k 
in Theorem 12 . 7. 2 1 the formula for averages of characteristic polynomials takes the following form: 



(2.8.4) 



det 



1 d(ui) 



Ui — Vj d(vj ) 



det 



m=V A N (f) 



1 



d(uj) 

Ui - vj \d(vj) 



'J > I A N (fl 



We multiply the right-hand and left-hand sides of the equation above by the product Yl ( u j ~ v j)- 
After that we differentiate k times with respect to u%, . . . u^, take u\ = x>\, . . . , Uk = Vk and find 
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(after simple manipulations) 

1 



(2.8.5) 



vi - v 2 
G(v 2 ) 



G(vi) 
1 

v 2 - Vl 

1 1 

Vk - Vl v k - v 2 

< G Oi)>A w (/) 

1 /d(v 2 )\ 
v 2 -vi \d(vi)/ 



V\ - v k 
v 2 - v k 



G(v k ) 



Aiv(/) 
1 /d(vi)\ 
v 1 -v 2 \ d(v 2 ) J 



Aiv(/) 



1 /d(v k )\ 
v k -v 1 \tZ(ui)/ 



A W (/) 



Ajv(/) 



where we have introduced 



1 /d(v k )\ 
v k -v 2 \d(v 2 )/ ANif) 

N 1 

g( v ) — X, v _ x . • 



3=1 



We note 



Res... Res (G(v i) . . . G(v k )) A (/) 

^1=3/1 V m =ym I " yJ > 



Xe Conf N {X) 



1 Mm)\ 

wi - v k \d(v k ) / 
v 2 - v k \d(v k ) J 



Aiv(/) 
Aiv(/) 



< g H)>a„(/) 



Prob(X) 



Qm(yi,---,Vm) 



Thus the residue of the left-hand side of equation Q2.8.5|) is precisely the m-point correlation 
function. Now we compute the residue of the right-hand side of equation (|2.H.5jl . It follows from 
the definition of K(x,y) (see equation (|2.8.3|) 1 that 



Res 

v=y 



(G(v)) 



Aat(/) 



K(y,y) 



The determinant in the right-hand side of equation (|2.8.5f) is a sum of products, and the off-diagonal 
elements are included as the cyclic products. The residues of the cyclic products are then the cyclic 
products of the kernels K(yi,yj), for example 

1 1 /d(vi)\ /d(v 2 )\ /d( Vj )\ 

yi-v 2 v 2 -v 3 vj-vi \d(v 2 ) / AN{f) \d(v 3 ) / AN[f) \d(wi)/ Ajv(/) _ 



Res . . . Res 

vi=yi Vj=yj 



K(yi,y 2 )K(y 2 ,y 3 ) . . . K(y j ,y 1 ) 



It means that the residues of the right-hand side of equation (|2.8.5j) is det [K(yi, J/j)]™- =1 which 
proves the Proposition. □ 



2.9. Proof of Proposition 12.5.31 This standard result (see Mehta Deift [22]) follows from 
the formula for the averages of two characteristic polynomials [341 1531 |5| , 



J_ /d(x)\ 
-t \d(0/ 



7iv-i 



Ajv(/) 



TT N (x)h N _ 1 (^) - TT N -l(x)h N (£) 
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Indeed, 



Res [Me)] = -7T - fivWiy) 



Therefore Proposition 12.8,51 gives 



when x / y, and 



K(x,y) 



K(x, x) 



f(y) 



Tr N (x)TT N _ 1 (y) -7rjv-i(x)7rjv(y) 



-JV-1 



x-y 



f{x) [7r^(x)7rAr_i(x) -7ij V _ 1 (a;)7rjv(a;)] 



"JV-l 



N-l 



/(*)£ 

j=0 



irj(x)7Tj(x) 



3 



when x = y. Taking into account the Christoffel-Darboux summation formula (see [HE] ) we see that 
the determinant with the kernel K(x, y) is equal to the determinant with the Christoffel-Darboux 
kernel K^ D (x,y). 

3. PFAFFIAN POINT ENSEMBLES 

3.1. Pfaffian L-ensembles. Given a set X let us construct two copies of X, and denote them 
by X' and X". We then introduce a 2 x 2 matrix valued function L(x,y), which depends on two 
variables x,y £ X 

~ L(x',y') L(x',y") 
L(x",y') L(x",y") 

Here we denote by the same letter L 2 x 2 matrix valued function of two arguments x,y £ X, 
and the scalar function whose arguments are taken from X' U X". Once x,y take values in X, 
the variables x' , y' (x" , y") are the elements of X' (X") corresponding to x, y. Assume that this 
function L is antisymmetric. Then L(x,y) defines a 2 x 2 block antisymmetric matrix on X. 

To any X C X there will correspond a 2 x 2 block antisymmetric submatrix of L. We denote 
this submatrix by L(X\X). If X consists of m points, 

X = (xi, . . .,x m ) , X C X 

the submatrix L{X\X) has the form 



(3.1.1) 



L(x,y) 



L(X\X) 



Lix'^x'l) 
-Lix' x ,x'{) ^ 



-L(xi,x m ) 
-Lix 1 ,x m ) 



-L{x'{, x' m ) 
-L(x± , x m ) 



L(xi,x' m ) 
L(x 1 , x m ) 



— Z/(x m , x rn ) 



L(x'[, x n 







Denote by Pf A the Pfaffian of an even dimensional antisymmetric matrix A. The definition of 
Pf A is given in Appendix, Section 16.61 Assume that the matrix L has the property 

PfL(X|X)>0, VXcX. 

Let J be 2 x 2 block matrix of format X x X with matrix elements 

1 

(3.1.2) J(x,y) 



1 



, x = y; 
0, otherwise 
Definition 3.1.1. A point process on X defined by 

(3.1.3) 

is called the Pfaffian L-ensemble 
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The fact that ^ Probi(X) = 1 follows from the expansion of Pf (J + L) into a sum of Pfaffians 

xcx 

of the symmetric 2x2 block submatrices L{X\X) of L, i. e. 

(3.1.4) Pf( J + L) = pf L{X\X). 

xcx 

The striking property of the Pfafhan L-ensembles is that the m-point correlation function g m (Y) 
is given by a Pfafhan, 

(3.1.5) Qm(Y)=Pf\K( yi , yj )]™ =1 , Y = ( yi ,...,y m ) 
Here the matrix K is defined in terms of L by the expression 

(3.1.6) K = J + (J + L)- 1 

the Pfaffian expression for m-point correlation functions reflects the fact that the Pfaffian L- 
ensembles is a special class of Pfaffian point processes. Different Pfaffian processes were considered 
previously by Tracy and Widom Rains Olshanski Soshnikov [5T] . 

3.2. Special matrices L. Here we assume that sets X, X are ordered. Thus if X = (xi,X2, ■ ■ .) 
then x\ < %2 < .... For any x E X we denote by l x the immediate left neighbor of X and by r x 
the immediate right neighbor of x. For example, if X = Z then ^x = x — 1 and r x = x + 1. 

Given a fixed splitting of X into positive and negative parts, X = X_ U X + , we denote by y 
the minimal (left most) element of X+. We introduce the parity on the sets X±, referring to the 
minimal elements of these sets as odd elements. 

According to the decomposition of the set X , X = X_ U j: U X + \ r, we write the matrix L in the 
block form: 







' L__ 


L-o 


L-+ ' 


(3.2.1) 


L = 


Lo- 


Loo 


Lo+ 






. L +- 


L+o 





We are interested in the matrices L defined by 
(3.2.2) L = 



E A 
A T 
B T 



B 





As usual, here E, A, B are the matrices with 2x2 block elements. Specifically, 



(3.2.3) 



e(x,y) 




x, y E X_ 



(3.2.4) 



A(x,y) 



e(x,y) 

h(x)h(y) 



x-y 



x E X_, y = p. 



(3.2.5) 



B(x,y) 





h{x)h(y) h(x)h{ l y) 



, x E X_, y E X+ \p. 



x — y x — L y 

The two-point function e(x,y) in equations above is antisymmetric, e(x,y) = — e(y,x). When 
x < y, 



(3.2.6) 



e(x,y) 



1, x — odd, y — even 
0, otherwise. 



The function /i is nonnegative on X. 

Configurations X E X can be divided into two classes. The first class consists of configuration 
which do not include the point p. Such configurations have the form X = X- U X + , X + = 
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(x^ ,x 2 , ■ ■ ■), %i > P- The second class consists of configurations including the point p. For such 
configuration X + = (p, xf, x 2 , ■ • •)• For any X C X denote by X the configuration defined by 



X_ = X. 



X 4 



x+np = 
x+m^ 



//++/++ \ 

I x l > x l ) x 2 ' x 2 ' ' ' v> 
^, , J,j , x 2 , x 2 , • • ■), 

Definition 3.2.1. We say that X G Conf L (X) if 

• X+ = (x+ < x\ < . . .) 

• all points of X + are different 

• |X_| = |X+| 

• X_ = (x^ < x 2 < • • •), where x~ has the same parity as i. 
This definition is justified by the following statement. 

Theorem 3.2.2. With L given by equations \H.2.1]) - \H.2.b}) we have 

1 V(X„)V(X + ^ 



(3.2.7) 



Prob L {X) 



Pf(J + L) n(X+;X_) 



h(X) 



for X G Conf L (X) and for all other X £ X. 



Proof. The positive integer d = \X- \ = |X+| can be even or odd, depending on whether X includes 
the point p or not. According to that consider two cases. 
Case 1. X n p = 

Given copies X' , X" of I G X in I', X" we denote by X' tb) X" the set {x'^x'l, x' 2 , x' 2 ', ...). Then 
we have 



Pf L{X\X) = Pf L 



Pf L 



I_UI,|I_UI, 



XL a XL) u (x; a x'i) |(xL « xL) u (xL a xr 
X X' 



d(d— 1) 

(-)-^ -PfL 



x Pf L 



x"_ u (x+ u x") |xL u (x+ a x+; 



as the function L(x, y) = for any a; G XL and any y which does not belong to X'_ (see equations 
(fTPTl - lEmi) and We note that L(x,y) = 0, if x,y G XL, or if x,y G X^ W X^ 



Therefore XLI = XLI + 



Consider Pf L 



XL IX' 



XL I, or X-l = 2IX4-I, which means that IX+I = |X_|. 



. Note that the matrix L 



XL I XI 



is even dimensional, if |X_| = 2|X + |. 



Moreover the matrix L 



XL ix! 



is the matrix whose (i, j) entry is, by definition, given by e(x i ,Xa ) . 



Clearly, if x 1 is even, the first row of this matrix consists of only zeros. Thus, if Pf L 



X' IX' 



^0, 



x x must be odd. Now assume that x 2i _i and x 2i have the same parity. In this case (2i — l) s ^ and 



th 

2i rows of the matrix L 



XL IX' 



are equal to each other. Therefore, if Pf L 



X' |X' 



7^ the 



elements of the set X_ = (x x ,x 2 , . . .) are such that x 1 is odd, x 2 is even, x 3 is odd and so on. 
This proves the condition on the parity for the configurations in Conf L (X). Moreover, using the 
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definition of Pfaffian it is not hard to conclude that Pf L 
non-zero probabilities. 



X'_\X'_ 



1 for the configurations with 



Since \X'i\ = \X' + \ + \X'±\ the matrix L 
structure: 



XL u (x' + a x'l) \XL u (x' + a x'l) 



has the block 



Odxd Qdxd 
-Qdxd °dxd _ 



with 



Q 



dxd 



h{xi)h(xt) h(x^)h( l xt) Kxi)h(x+ 2 ) h{ Xl )h{ l x + d/2 ) 



X X 



Xi X-\ 



Xi — X 



d/2 



Xi — l x 



l„+. 



h(x d )h(x±) h(d l )h( l Xi^ 



d/2 



h ( x d )h{x+ /2 ) h(x d )h{ l x d/2 ) 



x d x \ 



x d x \ 



(d is even). Thus we have 



d(d-l) 

Pf L(X\X) = (-J-hr^Pf 



X d X d/2 



"Qdxd °dxd 



x d lx d/2 



det Qrfxd 



n(*+;*-) 

d(d-l) . d 

where we have used the formula for the Cauchy determinant, see Appendix. Noting that (— ) 2 +2 

d 2 

(— )~2~ = 1 (as d is even) we obtain the formula stated in the Theorem. 
Case 2. 

The proof is very similar. We observe that any configuration X has a form 



Then 



Pf L(X|X) = Pf L 



X = X- U r U X 4 



X_ U r U X + 1 X_ U y U X 4 



Pf L 



(x'_ ax")u (y', r") u (x' + w x^) | (xi w x^) u (y', y") u (x; a x^) 



d(d-l) 

= (-)"^ - Pf L 



X^,y'|Xi,y' 



x Pf L 



x^ u y" u (x; a x£) |x^ u f u (x; a x^) 



Clearly, \X'L \ = |?' U (X^ W X^) |, otherwise Pf L(X\X) = 0. Thus X_ consists of odd number of 
elements, and |X_| = |X + | = d, d is odd, and we repeat the same computations as in the previous 
case. □ 



3.3. Pfaffian L -ensembles. Given X = X~ UX+ consider a different splitting of X, X = £_ Ul + . 
Here X + = X + U £0 , = X- \ Xo , and the set Xo is a subset of £_ . Assume that Xo consists of 
25 rightmost points of X- . 

Define a new matrix L on X in such a way that it has the same structure with respect to 
the decomposition X = X- U y U X + \ y as the matrix L with respect to the decomposition X = 
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£_ Uf U J6 + \ p. (Here p stands for the minimal (left most) element of Then L is given by 



(3.3.1) 



L 



E 


A 


B 








-A T 














-B T 









Here E, A, B are the matrices with 2x2 block elements. Specifically, 



(3.3.2) 



E(x,y) 



e[x,y) 




, x,y G X. 



(3.3.3) 



A(x,y) 



e(x,y) 0^ 

h(x)h(y) 







x-y 



x € y = p. 



(3.3.4) 



B(x,y) 





h(x)h{y) h{x)h( l y) 



x-y 



x — l y 



, x G y G £+ \p. 



If /i is nonnegative on X the matrix L defines a Pfaffian ensemble on X which we call L-ensemble. 

We introduce the set of point configurations Conf L (X) in the same way as the set Conf L (X) 
for the Pfaffian L-ensemble was introduced. Namely, for any Z G X denote by Z the configuration 
defined by 



z_ = z_ 



Here Z± = Z n £± . 



(£, l z-i ,zi , l z 2 ,z 2 , ■ ■ ■), z + np/ 



Definition 3.3.1. We say that Z G Conf L (X) if 

• all points of Z + are different 

• |Z_| = |Z+| 

• Z_ = (z^ < < . . .), where z~ has the same parity as i. 
Theorem 13.2.21 savs that with L given by equations (|3,3.1j) - H3.3.4j) we obtain 

(3.3.5) Prob.fZ) = ™ h(Z) 

Pf(J + L) n(z+;Z-) 

for Z G Conf L (X) and for all other Z C X. Given Z G Conf L (X) we build from Z a con- 
figuration X by the particle-hole involution on Xo (see Section 12.11 for the definition) . We note 
that the corresponding configuration X does not belong to Conf L (X), since the configuration X 
is unbalanced with respect to the splitting X = X- U X + . Moreover, 

• \X_ \ - \X+\ = 2S 

• X + = ( l Xi,xf, l X2,x%,...), A+np = 

• X + = (];, l x± ,xf, l x 2 ,x^, ...), X + Hp / 

• A_ = (x^jX^" , • • •), x~ has the same parity as i, 
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and all points of X± are different. 

Define the weight h of the .L-ensemble in terms of the weight h of the L-ensemble by the formula: 



(3.3.6) 



h(z) 



Hz) n \z- 

ye Xo 

1 



y\ 



z 6 X_ \ Xq , 



o, 



n 



ye % 



y| 



z6L 



Note that this formula is identical to ()2.3.3|) except for the absolute values. We did not need 
absolute values in ()2.3.3|) because there the formulas only contain h 2 and not h. 



Proposition 3.3.2. With h(z) defined by equation kS. S.b}) and X constructed from Z by the 
particle-hole involution on Xq, 



(3.3.7) 



Prob;(Z) 



1 



h(X) \V{X-)\\V(X+)\ 



Pf (J + L) \V(Xa)\h(Xo) \[{X + ,X 



Z G Conf L {X) 



Proof. The probability of the configuration Z can be rewritten in terms of X by the same method 
as it was done for the case of determinantal ensembles, see Proposition 12.3. II □ 

3.4. Correlation functions of Pfaffian L-ensembles and averages of characteristic poly- 
nomials. The goal of this Section is to express the kernel K of the m-point correlation function 
g m (Y) in terms of averages of characteristic polynomials associated with Pfaffian L-ensembles. 
Namely we want to prove the Pfaffian analog of Proposition 12.4.21 

Introduce nonintersecting sets of complex numbers with nonequal elements, 

(3.4.1) a+ = (af,...,a£), a_ = (a7/ , . . . , a~) . 
Assume that ^ fl I = 0, and 

(3.4.2) k-m= \a+\ - |a_| = 2S, S £ Z. 

Similarly to the case of the determinantal L-ensembles we extend the definition of the matrices 
K, L as follows. Let a', a" denote two copies of the set a = a_ U a+. We add to L rows and 
columns parameterized by a', a", and then define new matrix elements of L in accordance with 
equations ()3.2.2|) - (|3.2.5I) . where we assume that a_ is added to and a + is added to X + . Then 
we consider the matrix 

(3.4.3) K(a U £|a U X) = J + (J + L(a U X|a U X))" 1 . 

Proposition 3.4.1. The Pfaffian of the symmetric (k + m) x (k + m) submatrix K(a'\a') of the 
matrix K\ a' ttl a" U X' ttl X" \ a' l±) a" U X' t±J X" ) can be given as a normalized average of the functions 



E(-,-) defined by equations \2.4-^\ ) with respect to the Pfaffian L-ensemble. Namely, 



Pf [K(a'\a')] = h(a) 



V(a-)V(a + ) 



(3.4.4) 



n(a_;a+) 

n("+;*o) 



Pf (J + L) 
Pf (J + L) 



h(X )\V(X ) 



E(a + ,Z) \ 
E(a-,Z)/ L 



Proof. Combining equations (|6.8.3|) and (|6.9,1[) of Appendix we obtain the following formula: 

1 



(3.4.5) 



Pf K[a'\a] 



Pf( J + L)£ Pfi 



xe x 



a" UX' \t)X"\a" UX' ttl X" 
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Here a is a set with even \a\ = k + m, and a' , a" are two copies of a. Now we compute the Pfaffian 
in the sum using similar arguments as in the proof of Theorem 13.2.21 

Assume first that X n y = 0. We denote L(/3) = L{[3\(3) for any set j3. Then the Pfaffian in the 



sum is 



Pf L 



a" Ul'HJ X" 



Pf L 
Pf L 

(") 



a'L u a\ nx'_^x"_nx' + ^x'{ 
a'L uxiwx!uii.yx!u a" 

X_|(|X_|-1) 



■(_)l<*-ll*-lpf L 



Xi U a'L U X^ U Xl tt) X" U a'\ 



(-)■ 



|X_|(|X_|-1) 



>|a_||X_|p£ ^ 



Pf L 



a" uz!uxityx!u a': 



By equations (|3.2.2|) - (|3.2.5|) *+) = 0. Thus, the last expression does not equal to zero only if 

| a— | + |X_| = |a+| + 2|X+|, and X_ = (x7 ,x^ , ■ ■ •) is such that x7 has the same parity as i. Since 
|a+| — |a_| = 25", X_ consists of even number of elements. We note that Pf L{X'_) = 1. Then 



Pf L 



a" UX'W X" 



|X_|(|X_|-1) 



IJf+l 



V(a_ UX_)F(X+ Uaq 
\[{a UX_;X+ Ua + ; 
x h(X)h{a) 



(-)" 



|X_|(|X_|-1) 



'(-)- 

y(a_)y(a + ) 

]!(«-;«+) 



£(a_,X) 



y(x„)y(x 4 



h(X) 



It remains to see that the sign cancels out. We know that |X_| — |X_j 
|X_| is even. Then 



|a+| — |a_ 



25 and 



and the result is 



|X_|(|X_|-1) t |X + | 

= {-f 



(-) 

V(X„)V(X, 

U(x + ;x ; 



■(-)■ 



5 |V-(X_)||V-(X + )| 



Pf L 



q" U X'ttX" 



V(a-)V(a+) 
Yl(a-;a+) 



E(a + ,X) 
E(a-,X) 



\V(X_)\\V(X, 



MX). 



Now we apply Proposition 13.3.21 and obtain formula l|3.4.4|) . The case X n y ^ is considered in 
the same way. □ 
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3.5. Discrete symplectic ensemble. Given X denote by Conf^iX) the following set of point 
configurations: 

Conf^(X) = S.X C X\X = ( l Xl < x x < . . . < l x N < x N )\ 

Assume that a nonnegative function / is given on X, which does not vanish at least at 2N distinct 
points. 

Definition 3.5.1. The point process which lives on Conf^{X) and for which the probability of 
a configuration X is given by 

-l N 

(3.5.1) Prob(X) = [c$] Hf( Xi )\V(X)\, X e Conf^(X) 

i=l 

will be called 2N point discrete symplectic ensemble and will be denoted by A^(/). Here cff is 
a normalization constant. 



Remark 3.5.2. Rewrite |V(X)| as follows 



(3.5.2) \V(X)\ 



N 



In the continuous limit the points l x and x get closer and closer to each other, and all differences 
l X{ — Xi turns into the same small constant, say e. Thus the degree of the Vandermonde determinant 
turns into four, and the probability distribution 1)3.5.1(1 takes the same form as the probability 
distribution of eigenvalues for the symplectic ensemble of the Random Matrix Theory (see, for 
example, Mehta (44] . Chapter 3). 

For a symmetric function g(X) = g(x\, x\, . . . / xjy, xn) of points of the configuration X G 
Conf^(X), the average with respect to the discrete symplectic ensemble is defined by 

(3-5.3) <0) aW = [c^r 1 Y, g(X)\V(X)\f(X) 



Xe Conf^(X) 



If X is a lattice then in the continuous limit we introduce the skew symmetric inner product for 
arbitrary functions gi, g2 on X 

(3-5.4) ( gi ,g 2 ) = £ ( gi { l x)g 2 (x) - g 2 { l x) gi (x)\ f{ l x)f(x) 

where X denotes X without the leftmost point. If 

(3-5.5) det 11(^)11^/0 

then a family of monic skew orthogonal polynomials associated with the discrete symplectic en- 
semble can be constructed. 

Definition 3.5.3. For i = 0,1,2,... let p 2 i, P2i+i be monic polynomials of the degrees 2i and 
2i + 1, which satisfy the conditions 

• (P2i,P2j) = 0, (p 2i+ l,p 2 j + l) = 

• (P2i,P2j+l) = hiSij. 

The family {p 2 i,p 2 i+i} will be called the family of the skew orthogonal polynomials. 
Remark 3.5.4. The skew orthogonal polynomials are defined up to the replacement 

P2j+i P2j+i + const p 2j . 
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Lemma 3.5.5. For arbitrary functions (j>i(x), i = 1,...,2N and an antisymmetric two point 
function e(x,y) the following identity is valid 



(3.5.6) 

where 
(3.5.7) 



E pf 

xi,...,x 2 n£ £ 





2N 




2Af 




2N 


6 (Xi , Xj ) 


det 




= (2JV)! Pf 



















J l1 Yj/ 



^2 e(v,y)<i>i( x )<l>j{y) 

x,yG X 

This is a well known de Bruijn identity, see e.g. de Bruijn |19j . Tracy and Widom [HE], Baik and 
Rains _7 , Rains |49| . This statement readily implies, (see Rains |49| . Corollary 1.3): 

Corollary 3.5.6. Let (pi, ... , 4>2N ; ^l, • • • ; V*2JV & e arbitrary functions on X. We set 



and introduce the matrix A 
(3.5.8) 

Then 
(3.5.9) 







" M-) 








) = 




, #) = 








_ 027v(O 




_ ^2JV(-) _ 






2N 






























= E M^jfa) - 4>j 





x& X 



E det 



Xi<...<a;jy 
XjG X, i=l,iV 



Pf A 



The following statement is well known as well, and we give a proof for the reader's convenience. 

Proposition 3.5.7. The normalization constant cff is equal to the product of {p2i,P2i+i), i>e. 

JV-l 

(3.5.10) c% ] = Y[ h = (p 2 i,P2i+i) 



8=0 



Proof. Set 



TTi-l(x) = X 1 + . . . , i = 1, 



,2iV. 



This gives a system of polynomials of degrees 0, . . . , 2iV — 1 with the highest coefficients equal to 
one. The constant is the average over point configurations which includes the absolute value 
of the Vandermonde determinant and the product of weights. Rewrite the absolute value of the 
Vandermonde determinant in terms of the polynomials irj. In particular these polynomials can 
be chosen to be the skew symmetric orthogonal polynomials with respect to the weight / in the 
average. Then we obtain 



(3.5.11) 



,(4) 



E det 



X1<...<XJ£ 

Xi e X, i=l,N 



^(xi),ip(xi),... ,<P(xn),i/j(x n ) 



where we have introduced 

(3.5.12) <k{x) = Pl - 1 ( l x)f( l x), 4i(x)= Pi - 1 (x)f(x), i = TJN 
Now we apply Corollary 13.5.61 and prove the Proposition. 

Given X E Conf^(3i) and a complex parameter £ we define the characteristic polynomial 

(3.5.13) d(C) = n(C;*)- 

The following analog of Heine's identity is also well known, see e. g. Eynard Forrester |29j . 



□ 
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Proposition 3.5.8. 

(3.5.14) 



P2N(() = <d(C)> A (4) 



Proof. The average is a sum over the variables xi, . . . ,xn which take values in X. Set 

(3.5.15) <f>i(x) = Pi - 1 { l x)f{ l x), ip i (x)=p i - 1 (x)f{x), i = l,2N + l 

Then the expression under the sum is equal to the determinant of size 2N + 1 x 2N + 1, whose odd 
columns are 4>(xi),i = 1,N, the even columns are if)(xi),i = 1,N, and the last 2N + 1 th column 
is Pi-i(C)- We represent this determinant as a sum over permutations a € SW+i- Note that the 
sum over all permutations which do not satisfy the following two conditions 
• cr(2iV + l) = 2N + 1 

. | {(7(1), a(2)}, . . . , {a(2N - 1), a(2iV)}} = {{1, 2}, . . . , {2iV - 1, 2iV} 

vanishes after averaging. This is so because of the orthogonality relations for the skew symmetric 
orthogonal polynomials. On the other hand, there are N\ permutations which do satisfy the con- 
ditions above. All of them give the same contribution, which is proportional to P2n(()- Obviously, 
(d(Q) is a monic polynomial. This completes the proof. □ 

Proposition 3.5.9. For £ X the average o/l/d(£) is given by the skew symmetric inner product 
ofR^x) = ((-x)' 1 and p 2 N-2 (x), 



(3.5.16) 



<iMC)> 



(/) 



h N \{p2N-2,Rd- 



Proof. The average over the discrete symplectic ensemble is the normalized sum over configurations 
from the set Conf^l {X) . The normalization constant is cff . This normalized sum is that over the 
ordered variables X{ G X, i = 1,N. We have a symmetric function of variables xi under the sum, 
so we can remove the ordering of the variables changing the normalization constant from cff to 
N\$. 

The key observation which helps to compute the average is the following one. The expression 
under the sum contains the absolute value of the Vandermonde determinant, the product of weights 
and l/d((). This expression can be simplified if we expand 1/g?(C) into partial fractions, 

N 



^(0 = E 



(3.5.17) 



u=l s 



N N 

n ( i x u - Xj) n { i xv - i xj 

j=l 3=1 



N N 

n (xu -xj)u [x v 

3=1 3=1 
3^v 



_ I 



Xj) 



We note that each term of the first sum and each term of the second sum gives the same contribution 
to the average. Indeed, we always can make the change of variables under which x u becomes x±, 
and l x u becomes l x\, thanks to the symmetry of the involved functions under permutations of 
variables. Therefore the average will consist of two terms. The first term is 

£ ~N Rd l xi)d( Xl ) 



Xl,...,X N 



V( l x 2 ,x 2 , ■■■! x N ,x N ) 



N 



where we have denoted d(x\) = \\ {x\ — Xj)(x\ — l Xj). The second term is 

J'=2 



E N R c ( Xl ) d( l Xl ) 



Xl,...,X N 



V( l x 2 ,x 2 , • • • / x N ,x N ) 



f{X) 
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The computation of the sums is reduced to the computation of averages of characteristic polyno- 
mials over configurations X £ Conf^_ 2 {3i)- The previous Proposition says that these averages 
are skew symplectic polynomials of the degree 2N — 2. The remaining sum over x\ gives us the 
skew symmetric product in the righthand side of formula (|3.5.1(ij) . □ 

Given the family of the skew symmetric orthogonal polynomials with respect to the weight / 
introduce the Christoffel-Darboux type kernel 

1 P2i+1 (C)P2i (V) ~ P2i+1 (v)P2i (C) 



(3.5.18) K^%r)) = J2 



hi 

i=0 



Proposition 3.5.10. For £,7/ X, the averages of products and ratios of two characteristic 
polynomials with respect to the discrete symplectic ensembles are given in terms of K N D,4: and its 
pairings with R^ = (( — x) , R v = (?? — x) : 

(3.5.19) (d(0d(r))) A m w = ^ KS+?U,V) 

(3 ' 5 - 20) (ty). S( /) = {dMd{ M^ • m) 

) < «> AS(/) - ^ (*<•> (M)^_ M ^ 

Remark 3.5.11. Clearly, the first equation remains valid if C,7] £ X, and the second equation 
remains valid for rj £ X. 

Proof. Basically we proceed as in the computations of averages of the characteristic polynomial and 
its inverse. Let us first prove the formula for the average of product of characteristic polynomials. 
This average can be represented as a sum over the variables Xi £ X,i £ 1, iV, divided by the 
normalization constant N\. (The ordering of x±,X2, . . . ,xn is removed by AH). Set 

4>i(x) = P i„ l ( l x)f( l x), ^ < (x)=jii_i(x)/(x) ) i = l,2N + 2 

The expression under the sum can be rewritten as the determinant of the 2N + 2 x 2N + 2 matrix 
whose first column is Pi(C), the second column is Pi{rf), all other odd columns are (fi(x2i-i), and 
all other even columns are ip(x2i)- The average now is equal to this sum multiplied by the factor 

Now we rewrite the determinant as a sum over permutations a £ S2N+2, and change the order 
of sums. The skew symmetric product of (f>2%-\ and ip2i gives foj-i, and all other skew symmetric 
products constructed with these functions are zero. Then the sum over all permutations a which 
do not satisfy the relation 



j{a(3), a(4)}, . . . ,{a(2iV + 1), <r(2iV + 2)} J 



{l,2},...,{2i-l,2i},...,{2iV + l,2iV} 

for some 1 < i < N vanishes after the averaging. Consider the permutations which do satisfy this 
condition. All these permutations transfer the pair (1,2) to the pair (2i — l,2i). The sum over 
such permutations converts the sum over x%, . . . ,xn into the factor AH cisf+i/hi, thanks to the 
orthogonality relations between </> and ip. The index i takes values from 1 to A" + 1, and we obtain 
equation (|3.5.19j) . 

The second and the third equations are obtained from the first one by the same procedure as 
that used in the computation of (l/<i(£)) (4) m , Proposition 13.5.91 □ 
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Definition 3.5.12. The Cauchy type transform hk(Q of the monic skew orthogonal polynomial 
Pk is the skew symmetric inner product of and pk, 

(3.5.22) h k (C) = (pk,Rc) 

Here £ G" X, and the skew symmetric product is defined by equation (|3.5,4j) . 

Proposition 3.5.13. The averages of two characteristic polynomials with respect to the discrete 
symplectic ensemble are Christoffel-Darboux type kernels constructed from the skew orthogonal 
polynomials associated with this ensemble, and their Cauchy type transforms: 

• (d{C)d(v)) A w (f) = t — - E t 

/d(r,)\ , /~\ N ^ P2i + Mh 2i {Q-pK{rj)h 2i+1 {Q 

•W/^(/) = (7? "° S K + 1 



l \ 11 



\d(n)d(() / A w {f) h N _in-Q 



h^^h^C) - h 2t ( V )h2i + i(Q , , D p \ 

<L 7 1" \Rr]i K ti 

. i=0 n « 



Proof. Straightforward computation of the skew symmetric products in the righthand sides of 
equations ^EM ' ^M - □ 

3.6. Discrete orthogonal ensemble. Any point configuration from the set Conf^{X) defines 
a configuration of holes on X. All such configurations of holes naturally form a new set of point 
configurations. We will denote this set by Conf 2 ^{X). Thus, 

Conf^>(X) = |y CX\Y = X\X,X G Conf^(X) 

Point configurations from Conf^ admit an independent description. Namely, recall that X is the 
ordered set, in which the smallest point is the leftmost one. If we say that the leftmost point of the 
set X is odd, any configuration from the set Conf!^{X) is such that its smallest point is odd, and 
any two neighboring points are always have different parity. Assume once again that a nonnegative 
function / is given on X, which does not vanish at least at 2N points. 

Definition 3.6.1. The point process which lives on Conf^iX), and for which the probability of 
a configuration Y is given by 

(3.6.1) Prob(y) = const f(Y) ■ \V(Y)\, Y G Conf^(X) 

will be called 2N point discrete orthogonal ensemble and will be denoted by A^(/). 

Define a two point antisymmetric function e(x, y) on X such that it can be equal only to one, 
zero, or minus one: for x < y 

1, x is odd; 
0, otherwise. 



(3.6.2) e(x,y) 



Then the definition of A^(/) can be rewritten as 



(3-6.3) Prob(y) = [cf\ * f{Y) ■ \V{Y)\ Pf [e(y i} Vj )]» =1 

where Y is any configuration on X which consists of 2N points. This is the consequence of the 

r(l), 



-1 



fact that the Pfaffian in equation (|3.6.3j) is zero for any Y Conf^{X), and it is one for any 
Y G Conf^X). 

Remark 3.6.2. In the continuous limit the discrete orthogonal ensemble turns into the orthogonal 
ensemble of the Random Matrix Theory (see Mehta [44] . Chapter 3, for the definition). This is so 
because the parity condition for the configurations with nonzero probabilities becomes irrelevant 
in the continuous limit. 
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For any functions 51,52 introduce the skew symmetric inner product: 

(3.6.4) {91,92}= ^2 e ( x ii x 2)9i{xi)g2(x2)f(x 1 )f(x 2 ) 

£1,226 x 

and the monic skew symmetric orthogonal polynomials pi with respect to this inner product. (We 
assume that the weight / is such that det||(a; l ,a; j! )||fj_ 1 7^ 0). These skew symmetric orthogonal 
polynomials are constructed so that they satisfy the same orthogonality conditions as in Definition 
13.5.31 but with respect to the skew symmetric inner product defined by equation (|3,6.4j) . Word- 
for-word repetition of arguments of the previous Section gives 

Proposition 3.6.3. The normalization constant and the averages of characteristic polynomials 
must be of the same form for the discrete orthogonal and symplectic ensembles, up to the definition 
of the skew symmetric inner product. In particular we have 

N-l 

(i) 



(3.6.5) 4/ = Y\ h i, hi = (p2i,P2i+i) 



i=0 



(3-6.6) p 2N (C) = WC)) A (i) (/) 



(3-6.7) ( d_1 (C)) A W (/) = h N \{p 2N ^R c ) 

(3-6.8) < d (0«ffa)> A W w = ^ *$+}{C,V) 



/ 1 \ /.mW 



Here the Christoffel-Darboux kernel K NJ \ (C v) ^ s given by formula h3.5.18^ . where the skew sym- 
metric orthogonal polynomials pi are constructed with respect to inner product \3.6.4\l - 

Define the Cauchy type transforms of the skew orthogonal polynomials by the same way as they 
were defined in the case of the discrete symplectic ensembles, see Definition 13. 5 . 1 2l Then we obtain 
the following 

Proposition 3.6.4. The averages of two characteristic polynomials with respect to the discrete 
orthogonal ensemble are Christoffel-Darboux type kernels constructed from the skew orthogonal 
polynomials associated with this ensemble, and their Cauchy type transforms: 

IMS\M \\ Hn P2i+l(C)P2i(»7) ~ P2i(C)P2i+l(v) 

• (d(Qd(r])) w = - — - 7 

^nU) c, — J) i=0 hi 

\m/^ n = ( " " £ * + 1 



KM) 
(->—) 

\ d(v)d{Q I A m (f) h N -i 77 - C 



•y Wfo)MC)-Mi7)h2i+i(C) 1 m b\ 

2_. 7 1" {Mv,Mo 

. i=0 n i 



Proof. Straightforward computation of the skew symmetric products in the righthand sides of 
equations ilTTTs!) - itTTTTTiT. . □ 
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3.7. Relation of ensembles. This Section discusses the relation between Pfaffian L-ensembles, 
discrete symplectic ensembles and the discrete orthogonal ensembles, introduced in Section 13.21 
Section 13.51 and Section 13.61 respectively. Given X consider the Pfaffian L-ensemble, defined by 
equations (JS^Tjl - ipHO^l . Denote \X-\ = 2M, \X + \ = 2N, M,N £ Z >0 , and assume that the 
weight h is strictly positive on X. Given h introduce the weights by 



(3.7.1) 



/ {4) (^ 



h(x) 



x £ X^ 
x G X- 



(3.7.2) 



/ (1) (x) 



Uyex+ \ x ~_y\ 



, M^) n^/e X+, y+x \ x - V\ ' 



X £ X- 
x £ Xa 



Note that the only difference in the definition of /W and /( 4 ) is the interchange of X- and 3C+. It 
follows from the two equations above that 



(3.7.3) 



(4), 



n 

y& X, y^x 



X £ X 



Theorem 3.7.1. The Pfaffian L-ensemble with the weight h defined by equations IH. 2.1\) - i3.2.b}) . 

the discrete symplectic ensemble with the weight f^ defined by equation \H.5.1\) . and the discrete 
orthogonal ensemble with the weight f^ 1 ', defined by equation 113.6.1)) are equivalent in the follow- 
ing sense: 

a) There is a bijection between the sets Conf L (X) and Conf^{X) defined by the particle-hole 
involution on X+. Under this bijection the probability distribution for the Pfaffian L-ensemble, 
equation \3.2. 7| ), turns into the probability distribution for the discrete symplectic ensemble, equa- 
tion fTTHl) . 

b) There is a bijection between the sets Conf L (X) and Conf^{X) defined by the particle-hole 
involution on X_. Under this bijection the probability distribution for the Pfaffian L-ensemble, 
equation \3.2.T\) , turns into the probability distribution for the discrete orthogonal ensemble, equa- 
tion 13X71) . 

c) There is a bijection between the sets Conf^(X) and Conf^{X) defined by the particle-hole 
involution on the whole set X. Under this bijection the probability distribution for the discrete 
orthogonal ensemble, equation h3.6.1\) turns into the probability distribution for the discrete sym- 
plectic ensemble, K3.5.1]) . 

In other words, the probability spaces defined by Conf L {X), Conf^(X), Conf^{X) with the 
associated probability distributions are isomorphic. 

Proof. Given Conf L (X), let X £ Conf^iX), and X be the point configuration associated with X, 
see Section l3~2l Consider the particle-hole involution on X- . Under this involution the configuration 
X turns into the configuration Y related with the configuration X by the formula: 



(3.7.4) 



y = Y_uy+, y_ = x_\x_, y + 



x 4 



The key observation is that the configuration Y is in Conf^(X). It can be immediately seen 
from the definition of Conf^iX). Clearly, the particle-hole involution on X- is a bijection be- 
tween Conf L (X) and Conf^j^X). In order to prove a), it remains to show that the probability 

distribution on Conf L (X), equation (|3.2.7|) . turns into the probability distribution on Conf^ f (X), 
equation ()3.5.1|) . provided the weights are related by formula (j3.7,l|) . This is a consequence of two 
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facts. The first one is |V(Y")| can be rewritten as 

n n \x- y \ 

(3.7.5) \V(Y)\ = const V{X t )V{X - ) X+ V **~ , f , const = \V(X-)\ 

n(*+;*-) n n \%-y\ 

This follows from the relation Y and X (equation (|3.7.4|) ). The second fact is the relation: 

f(Y) = const ^ X+ \ const = /(£_) 

We then see that the righthand side of equation (|3.2.7[) is identical to that of equation (J3.5.1|) , if h 
and / are related by equation ()3.7.1|) . Thus we have proved a). 

The proof of b) is constructed in the same way. The configuration Y obtained by the particle-hole 
involution on X_ from X is 

(3.7.6) Y = Y- U Y+, K_ = X- , Y + = X+\X+. 
Instead of formula (|3.7,5|) we obtain 

(3.7.7) \v<r)\=™* V <*+> v f-> „ n^' 1 , - r . const _ \V( X+ )\ 

n(^+-,^-) n n \x-v\ 

which leads (together with the relation of the weights h and f^ l \ equation (|3.7.2j) ) to the equiva- 
lence of the probability distributions on Conf L (X) and Conf^iX). 

The statement c) obviously follows from a) and b). However, it is easy to give an independent 

argument. Let X £ Conf^^X). Then the particle-hole involution on X gives a configuration Y 
which is in Conf^iX). For X and Y related by the particle- hole involution we obtain 

(3.7.8) \V(X)\ = const \V(Y)\ \ [ f[ \ x ~Vl const = (y(X)r 1 

This expression (together with the relation of weights f^ and f^\ equation (|3,7.3j) ) leads to the 
equivalence of the discrete symplectic and the discrete orthogonal ensembles. □ 

Assume that the Pfaffian L-ensemble is given, and denote by A 2M (/( 4 )) and A 2 jv(/ (1) ) the 
equivalent discrete symplectic and orthogonal ensembles obtained from the Pfaffian L-ensemble 
by the particle-hole involution. With the Pfaffian L-ensemble we can construct a new ensemble, 
Pfaffian L-ensemble, by new splitting of X, X = X_ U X+, \X + \ = 2N + 2S, = 2M - 25, 
see Section 13.31 This Pfaffian L-ensemble induces a discrete symplectic ensemble of 2(M — S) 
point configurations, and a discrete orthogonal ensemble of 2(N + 5) point configurations. The 
first one is obtained by the particle-hole involution on X_, and the second one is obtained by the 
particle- hole involution on X+. Let h denotes the weight of the Pfaffian L-ensemble. Recall that 
it is given in terms of the weight h of the original Pfaffian L-ensemble by formula (|3.3.6j) . The 
Theorem above says that the Pfaffian L ensemble with weight h, the discrete symplectic ensemble 
A 2 (a/_5')(/( 4 )), and the discrete orthogonal ensemble A 2 (jv+s)(/^) are equivalent provided the 
weights are related by formulae (|3.7.1j) - Q3.7.2|) . with X± — > X±. 

Proposition 3.7.2. The discrete symplectic ensembles A 2 m(/ ( - 4 ^) and A 2 (m-S) 

(/W) have the 

same weight: = f^. The discrete orthogonal ensembles A 2 jv(/W) and A 2 /v+2S 

(/«) have 

the same weight: /W = f^ 1 '. 

Proof. Application of the same arguments as in the case of the determinantal point ensembles, see 
Proposition 12.6.21 □ 
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3.8. Averages of characteristic polynomials: symplectic and orthogonal ensembles. Our 

goal now is to derive averages of product and ratios for symplectic and orthogonal ensembles. We 
will use the relation of these ensembles with the Pfaffian L-ensembles discussed in the previous 
Section, and Proposition 13.4. II We begin with symplectic ensembles. 

Given equation (|3.4.4j) let us rewrite its righthand side in terms of the discrete symplectic 
ensembles. We proceed as in the determinantal case, see Section l2~Tl 

Proposition 3.8.1. The constant in equation \3.4-4\) * s the ratio of the normalization constants 
for the discrete symplectic ensembles of 2N — 2S particles and of 2N particles, 

(3.8.1) £^±|j = |f 

Proof. There are two Pfafnan ensembles involved in Proposition 13.4. II The first one, the Pfafnan 
L-ensemble, is defined with respect to the splitting X = X_ U £+, and the second one, the Pfaffian 
L-ensemble, is defined with respect to the splitting X = 3t_ U X + . The expression [Pf(J + L)]~ 



Pf( J + L) 



is 



is the probability of the empty configuration for the first Pfaffian ensemble, and 
the probability of the empty configuration for the second Pfaffian ensemble. 

By the particle-hole involution on X- and X_ construct the equivalent discrete symplectic en- 
sembles. Then the empty configurations of the Pfaffian ensembles turn into certain configurations 
of these symplectic ensembles. It is not hard to see that these configurations are and X-. 
Equate the probabilities: 

' ' f(3L.)\V(3L.)\, _J_ = -i_ 



Pf(J + L) C (J) J " " v Pf(J + L) c g)_ 5 

These equations, together with the relation between weights, equation (|3.7,1|) . give formula (|3.8.1|) . 

□ 

Theorem 3.8.2. For any integer N > 1 take an integer S such that 2 — 2N < 2S < \X\ — 2 — 2N , 
complex numbers a = (a_,a + ), a_ = (a^, . . . a~ ), a + = (af , . . . , at) such that \a+\ — |a_| = 
k — m = 2S, in each set a± the numbers are pairwise distinct, and the set a does not intersect X. 
Then the average of products and ratios of characteristic polynomials with respect to the discrete 
symplectic ensemble A^(/) is given by the formula: 



where the kernel function wff is defined by 

. W§ \a+, a}) = hxls^af - a+) (d(a+)d(a+)) 4) 



W^\a\a) 



< 4) (a-,a+) = -< 4) (a 



W { ^\a i ,a j ) = h N+s (a 




,(/) 



A (4) (f) 



a (4) m 



Proof. Assume first that S > 0. We use equation (|3.4.4j) . Denote by A^_ 2S (f^) the discrete 

symplectic ensemble obtained by the particle-hole involution on X from the Pfaffian L-ensemble 
involved in equation (|3.4.4j) . We want to rewrite the righthand side of equation (|3,4.4|) in terms of 
A^_ 2S (f^). A^_ 2S (f( 4 }) and the Pfaffian L-ensemble involved in equation (|3,4.4|) are equiv- 
alent, the corresponding probability measures are equal to each other. We only need to express 



:-!S 
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the functions E{. : .) in terms of characteristic polynomials associated with A^_ 2S (/^). For any 
finite set a, a n X = 0, we find 

E(a,X) - l[(a;X-) = ]J(a-Y) = Hd(ai) 



where the configurations X € Conf^(X) and Y E Conf!^_ 2S (£) are related by the particle- 
hole involution. This immediately follows from the definition of E(.,.), equation (|2.4.2[) . This 
observation and Proposition 13. 8. II give 

^( 4 ) 

(3.8.3) Pf [K{a'\a')] 



C »- s h{a) V{a - )V{a+) 



J4) 



II («+;%- 
n(a~;X_ 



n£i d K r )/ a 



;(/) 



Considering particular cases corresponding to k = 0, m = 2; k = 1, m = 1, k = 2, m = determine 
the kernel function K in this equation, and obtain the formula stated in the Theorem. 
The case 5 < corresponds to the splitting 



X = X_ U X H 
and is considered in the same way. 



X_ UX, 



0: 



X+ \ Xq 



□ 



Applying the same arguments we prove the corresponding result for the discrete orthogonal 
ensembles. 

Theorem 3.8.3. For any integer N > 1 take an integer S such that 2iV + 2-2|X| < 25 < 2N-2, 
complex numbers a = (a_,aq_), a_ = (aj~,...a~) ; a + = (af , . . . , at) such that \a+\ — |a_| = 
k — m = 25, in each set a± the numbers are pairwise distinct, and the set a does not intersect X. 
Then the average of products and ratios of characteristic polynomials with respect to the discrete 
orthogonal ensemble A^(/) is given by the formula: 



(3.8.4) 



nti d(at)/ A V(f) 



n(a-;a+) 



C W K(a_)7(a+) 



Pf 



W^ 1} (a|a) 



where the kernel function is defined by 



. W£\a+,a+) = h 



W; 



(1) 



A' 



• w 



(1) 

AT 



(«r» a j) 




,(/) 



hN-S-l 



2JV-25+2 



3.9. Computation of the m-point correlation function from the averages of character- 
istic polynomials. The aim of this Section is to illustrate the convenience and the generality of 
the averages of the characteristic polynomials. Similary to the case of the polynomial ensembles 
we extract the standard m-point correlation function from this averages. 

Proposition 3.9.1. The m-point correlation function of the discrete orthogonal ensemble A^(f) 
is given by the formula: 



m(m — 1) 

(3.9.1) 6m { Zl , ...,z m ) = (-)— 2— Pf 



[-5(Zj,Zi)]™. =1 



where 



h N . 



(C-rj) {d(0d(v)) A v 



,(/) 
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7(C,7?) = Res 
y=n 



Res 

y=r] 

Res 

y=v 



1 /d(0\ 
C-y \d(y)/ A W if) 



d 



dx 



-- y \ d (v) I A. ( 



&</>■ 



Proo/. Set q_ = (an, . . . , x m ), a+ = (t/i , ■■■,y m , 
determinant to rewrite equation Q3.8.4|) as 



det 



d(xj) 

Vj d{vj) 



C = V 

in Theorem l3.8.31 Use the formula for the Cauchy 
W$\x,y\x,y) 



m(m — 1) 



Pf 



« =1 / A«(/) 



Multiply the left-hand side of this equation by the product (xi — y{) . . . (x m — y m ). Differentiate 
the obtained expression by x%, . . . x m and take parameters Xj to be equal to the parameters yi. 
We then obtain some function of y\, . . . ,y m . By taking residues at points y\ = Z\, . . . ,y m = z m 
from this function we get Q m (zi, . . . , z m ) as we have learned previously in the determinantal case, 
Proposition 12.8,11 

We want to understand what happens with the righthand side under the same transformations. 

m 

Decompose the Pfaffian into the sum over permutations. Now we multiply this sum by Y\ ( x i ~ 

i=i 

differentiate and substitute xi = yt, i = l,m. We observe that all terms of this sum which do 

not include W^\xi, yi) remain unchanged, we only replace X{ by yi. As for the elements like 

W^\xi,yi) they turn into the elements 

W y N '{xi,yi) -> — {-rr-c ) 

dxi \ d(yi)/ A m {f) 



Thus the obtained sum is equal to 



Pf 



D(yi,yj) 


m 


S(yi,yj) 


m -. 

M=l 


-S{yj,yi) 


m 

i,i=i 




m 

i,j=l - 



where 



D{yi,yj) 



llN-l 

I(yiiVj) = h N {y% - yj) 
l 



{yi-yj){d{y i )d{y j )) w 



»(/) 



S(yi,yj) 



d(yi)d(yj] 
I d{yi 



y% - 

d 

dx 



Vj \d{yj) 

d{x) 
d{yi) 



2JV+2^ 

i + j 



It remains to take the residues. Rewrite once again the Pfaffian as a sum over permutations. Tak 
g into account the definition of the Pfaffian it is not hard 

the two point functions D, I, and S turn into D, S and I. 



ing into account the definition of the Pfaffian it is not hard to see that after taking Res , . . . Res 

y\=z\ y m =z m 
□ 



Proposition 3.9.2. The correlation functions of the discrete symplectic ensemble are given by the 
formulas of Provosition 9. 1\ with averages over A^ 1 ) replaced by similar averages over A^. 



40 



A. BORODIN AND E. STRAHOV 



Corollary 3.9.3. The m-point correlation function for the discrete orthogonal ensemble A^(f) 
is given by the quaternion determinant Tdet of the 2x2 matrix valued kernel a(zi,Zj), i,j = l,m, 



(3.9.2) 

where 

(3.9.3) 



i ■ ■ ■ i ^m) 



a(zi,Zj) 



Tdet a(zi, zj)™j =1 

~S(zi,Zj) I(zi,Zj) 
-D(zi,Zj) -S(zj,Zi) 



the two-point functions D, I, S are given by 



D(zi,Zj) 



N ^ P2l+l(zi)P2l(zj) ~ P2l(zi)p2l+l(zj) 



1=0 



I(Zi,Zj) = f{Zi)f(Zj) 



N-l 

£ 

L 1=0 



hi 

q2i + i(zi)q 2 i(zj) - q2i(zi)q 2 i+i(zj) 
hi 



+ e(zi,Zj) 



N-l 



£ 

1=0 



• S(zi,Zj) = f(zj) 

and the functions qt{z) are defined by 
(3.9.4) 



P21+1 (zi)q 2 i (zj ) - p 2 i (zi)q 2 i+i (zj) 



hi 



Res h k {() = f(z)q k (z) 

C=z 

Remark 3.9.4. With the specific choice e(x, y) = e(x — y) on a one dimensional lattice, 

' -1/2, x>y 

(3.9.5) e{x,y) = i 0, x = y 

1/2, x<y 

the formula for the m-point correlation function is reduced to that of Nagao and Forrester [U3 . 

Proof. Compute the two point functions D, I, S in Proposition 13. 9 . II We use the formulas which 
express averages of two characteristic polynomials as the Christoffel-Darbyoux type kernels con- 
structed from the skew orthogonal polynomials and their Cauchy type transforms, Proposition 
13.6.41 Taking into account equation (|3.9,4j) we immediately obtain the expressions for D, I, S 
stated in the Corollary. Thus we have 

Q m (zi, . . . ,Z m ) = (-) 2 Pf 

_ [-S(zj,zi)]™ =1 [I{z t , Zj )}™ j=l 

It remains to rewrite the righthand side of this expression as the quaternion determinant of 2 x 2 
matrix valued kernel. This is achieved in two steps. The first step is to observe that 



(3.9.6) Pf [W(x 1} . . . , x m ; yi, . . . , y m \xi, . . . , x m ; y u . . . , y m ) 
where 

V)(Zi,Zj) = 

The second step is to apply the formula 



m(m—l) 

(-)-^Pf 



w(Zi,Zj) 



»J=1 



W(xi,Xj) W(xi,yj) 
W(yi,Xj) It '!//,.//,) 



Pi[ZQ] = Tdet Q, Z 



(J ... \ 
J ... 



where Q is a 2 x 2 block antisymmetric matrix, and J 



\ ... J ) 

1 
-1 



□ 



Remark 3.9.5. The statement similar to Corollary 13.9.31 holds for the discrete symplectic ensem- 
bles as well. The proof is just the same. 
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4. Continuous Limit 

4.1. Exact formula for unitary (/? = 2) ensembles. Let \x be a positive measure on R with 
finite moments: J R |x| n //(dx) < oo for any n > 0, such that 

C ^ ):= M / II \xi-Xjfti<tei)---vi{dx N )?0 



i<i<j<N 



for any iV > 1 and some (equivalently, all) /? > 0. Denote by the probability measure on R N 
given by 

p ( fi\dx 1 ,...,dx N ) = Y\ \xi - Xjf fi(dxi) (g> ■ ■ ■ (g> fi(dx N ). 

N*-C N i<i<:j<N 

For any symmetric function g : E> N — ► C set 

(#) A (« := J g(x 1 ,...,x N )p ( fi\dx 1 ,...,dx N ) 

R N 

provided that the integral converges. For £ e C set 

£>(C) = II(C-^)- 
i>l 

The number of factors in such products will always be clear from the context. Recall that for finite 
sets A = {a\, . . . , a p } and B = {b\, . . . , b q } we use the notation YK^'i B) to denote the product of 
pairwise differences 

H(A,B) = f[f[(a i -b j ), 

i=ij=i 

and V(A) stands for the Vandermonde determinant ni<i<j< P ( a « — a j) (note that we use an ordering 
of aj's here). 

Theorem 4.1.1. (Case (3 = 2) For any integers N > 1 andS>l — N, and finite sets of mutually 
distinct complex numbers 

oT = {a^ . . ,a^J, a + = {af, . . . , a~}, 
(3- ={^,...,(3~ 2 }, /3+ = {/3+ 

with \a~ \ — \a + \ = \@~\ — |/3 + | = S, such that 

a~na + = 0, /Tn/3 + = 0, a + nM = 0, (3 + n K = 0, 

one /ias 

' n" i i^K-)ngf 1 ^r) \ _ r _ 1 ^[(i«-i+i/»-D a +i/»-i-i«-i] 



(4.1.1) 



,n?=i%;)nji^ 



cgjg n(»';« + )n(r ; /3 + ) 

C-W F(a-)F(«+)y(/3-)y(/3+) Qet 



^ 2) (a-,/3 + |/?-,a + ) 
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is a matrix with rows parameterized by elements of a and (3 + , columns 



parameterized by elements of (3 and or , and with matrix elements 

C (2) 



N+S-l 




N+S 



N ' ^N+S 

r (2) / 1 \ 

° JV4- .<?-!- 1 / -L \ 



Proof. Since all the moments of fi are finite, and we are averaging functions of at most polynomial 
growth, it suffices to prove the theorem for the measures with compact support. Clearly, formula 
(|4.1.1|) remains intact under the operation of taking weak limits of measures /i having the same 
compact support. On the other hand, formula (|4,1.1|) was already proved in Theorem I2.7.2I for 
all measures \i supported by a (large enough) finite set of points. Since any compactly supported 
probability measure can be weakly approximated by finitely supported ones, the proof is complete. 

□ 



4.2. Exact formulae for orthogonal (3 = 1 and symplectic = 4 ensembles. 



Theorem 4.2.1. (Cases (3 = 1 and (3 = 4) (i) For any integers N > 1 and S > 1 — N , and finite 
sets of mutually distinct complex numbers 



a = {at,...,a k }, (3 = {ft, . . . ,(3 m }, 
0, one has 



(4.2.1) 

such that (3 n 1 
(4.2.2) 

where wffl is a skew- symmetric (k + m) x (k + m) matrix with rows and columns parameterized 
by elements of a and (3, and with matrix elements given by 



n"i^(A)/ A ( 1 ) C$ V(a)V((3) 



Pf 



k — m = 2S, 



W£\a,f3\a,(3) 



Wjp(aii,aj) 



W^((3 h (3 3 ) 



a 



(i) 

2^+25-2 



(at - atj) {D(ai)D(aj)) ( i) 



rj( L ) A 2AT+2S-2 

1 /D(oq) 



Oi- (3j \D(f3j)/ A m 

^2^+23 



a 



a) 

2^+25+2 



a 



(i) 

2N+2S 



i 



r»(ft)r»(ft)/ A (x) 



2JV+2S+2 



(ii) For any integers N > 1 and S > 1 — N , and finite sets of mutually distinct complex numbers 
14.2.1)) such that (3 D K = 0, one has 



(4.2.3) 



C 



(4) 

N+S 



Pf 



wj}\a,P\a,p) 
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where wff is a skew- symmetric matrix with rows and columns parameterized by elements of a and 
B, and with matrix elements given by 

c (4) 

Wjpfaaj) = ^f" 1 (a t - ay) (Z> 2 (a i )D 2 (a i )) w 

Cn+s iV+s - 1 



a(4) 



C (A) ' I \ 



y Af+5 1 v v •" ' A n+s+i 

Proof. Arguing as in the proof of Theorem 14,1.11 we see that it suffices to prove Q4.2.2|) and Q4.2.3|) 
for compactly supported measures fi with continuous density with respect to the Lebesgue measure: 

H(dx) = F(x)dx, F e C(R), supp(F) c [-A, A], A G R+. 

Let us consider the case B = 1 first. The idea is to reduce the statement to Theorem 13.8.31 Split 
the segment [—A, A] into M > 2N intervals of length 2A/M: 



(4) 



-A, A] = Ji U • • • U Im, h 



, 2(fc-l)A a . 2fcA 
M ' M 



and take a discrete set X with 2M points such that each of the intervals Ik contains exactly two 
points of X. One may take for example 

-V — ( A i _A_ _A i 3A . _ a i M. _J i 14 ■ - A _ 3A 4 A.\ 

X. — \ 2M , /i -I- 2M , /i t 2M ' ' 2A/ i • • • ' 2M ' ^ 2M J ■ 

Then for any choice of 2N intervals J^., j = 1, ... ,27V, with fci < < • • • < fc2iv, there exists 
a unique way to choose a point of X in each of these intervals so that the chosen 2N points of 
X form a point configuration from Conf^{X). Conversely, any element of Conf^{X) produces 
2N intervals It which contain the points of the configuration. Thus, for any continuous function 
g:[-A,Af N 4c, _^ 

(H) g(x 1 ,...,x 2N ) 

(x 1 < — <x 2 N)&Conf^(X) 

can be viewed as a Riemannian sum for the integral 

g(x%, . . . ,x 2 N)dxi ■ ■ ■ dx 2 N- 

-A<X!<---<x 2N <A 

Taking M to oo brings the diameter of the partition [—A, A] = I± U • • • U Im to zero, and formula 
Q4.2.2JI directly follows from Theorem 13.8.31 The proof of (i) is complete. 

Let consider the case B = 4 now. We want to derive Q4.2.3j) from Theorem 13.8.21 Take M > 2N 
and 

, {2k-l)A\ M , . 
Then for any continuous function g : [—A, A] — > C, 



/ 



£ g(x 1 ,...,x N )\V(X)\ = (^y 



X=( l x 1 <x 1 <---< l x N <x N )£Conf^(X) 

^2 g(x!,...,X N ) Y[ ( x i ~ x j) 2 ( x i~ x j ~ w) ( X i ~ X 3 + W) 

Xi£X l<i<j<N 
-A+^<x 1 <-<x N 

= J g(xi,...,x N ) Y\ (xi - xj) 4 dxi- ■ ■ dx N + O (M^ 1 ) . 

~A<x 1 <-<x N <A l<i<3<N 



AA 



A. BORODIN AND E. STRAHOV 



Thus, Theorem ESU implies lUTST) . □ 

Remark 4.2.2. The conditions imposed on the measure [i at the beginning of this section can 
be relaxed. Namely, instead of finiteness of all the moments it suffices to require the convergence 
of the averages involved in (|4.2,2j) and (|4.2.3|) . and instead of requiring all to be nonzero it 
also suffices to require that only of the constants involved in these formulas. This allows to have 
finite support. 

Remark 4.2.3. The statements of Theorems 14.1.11 and 14 . 2 . T1 admit an interpretation in terms of 
random matrices. Indeed, the measures , A^y, and A^ are radial parts (i.e., projections onto 
(different) eigenvalues) of the measures 

P[dH) = const • exp(-Q(H)) dH, 

where Q(x) is an even degree polynomial with positive highest coefficient, and H belongs to the 
linear space of N x N Hermitian, 2N x 2N real symmetric, or N x N quaternion real Hermitian 
matrices, respectively. This corresponds to our measure \x being equal fx(dx) = ex.p(—Q(x))dx for 
[3 = 1 and 2, and fi(dx) = exp(— 2Q(x))dx for j3 = 4. In the cases = 1,2 the characteristic 
polynomial of the random matrix H, det(£ — H), is exactly our product-function D(Q, while for 
j3 = 4 we have det(£ — H) = D 2 ((). Thus, in Theorems 1,2 we deal with average values of products 
and ratios of characteristic polynomials with respect to the measures on matrices introduced above. 
In the case Q(x) = constat 2 these matrix probability spaces are called Gaussian Unitary (f3 = 2), 
Orthogonal ((5 = 1), and Symplectic ((3 = 4) Ensembles (GUE, GOE, GSE, for short). 

Remark 4.2.4. Formulas (|4.1,lj) . (|4,2.2j) . (|4.2,3|) require that the total number of factors in the 
left-hand side is even (we count factors in both the numerator and the denominator). These 
formulas can be easily extended to the situation when the total number of factors is odd. In order 
to do that, one needs to take the even case and send one of the parameters af,f3^ in (|4.1. If) and 

cti, f3j in (|4.2.2|) . ()4.2.3|) to oo. Given that D(Q ~ £ N in A^, the limiting formulas are readily 
obtained. After the limit transition, some of the two-point averages in the matrix elements of 
wffi in the right-hand side turn into one-point averages of the form (D^)), (1/D(j)) in the case 
(3 = 1,2, or (£> 2 (7)), (1/D 2 ^)) in the case = 4. 

4.3. Correlation functions. Theorems 14. 1 . II and 14 . 2 . ll allow to compute the correlation functions 
of A^ for (3 = 1,2,4 and reproduce fundamental results of the Random Matrix Theory. 
Define the nth correlation measure of A^ by 

Pn,N( dxi > • • • ' dXn ^> = (]v~n)! / PN ) (dx 1 ,... ,dx N ). 

X„ + l,X„ + 2,...,X N 

For a function f(z) of a complex variable C, which is continuous in both half-planes Im^ > and 
Im£ < up to the real axis, we will denote by [/(C)]^ =:r , ifl, the difference of the limit values 
of /(C) as C x from bottom and from top divided by 2ni: 

\f(0]x = ^ ri (f(x-iO)-f(x + iO)). 

Theorem 4.3.1. Take any n > 1 and assume that near n points x%, . . . ,x n £ M. the measure n 
is absolutely continuous with respect to the Lebesgue measure, and its density their is uniformly 
Holder continuous. Then the nth correlation measure of (or A^ for (3 = I) has a continuous 
density near (x\, . . . , x n ) which is given by 
• For (3 = 2 

pi 2 ]v(^i> • ■■,x n ) = det [K^\xi,Xj)]l j=1 
where for x ^ y the kernel is given by 



K^(x,y) = — 



D{x) \ 
D(C)/ 



A (2) 
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and for x = y the kernel is defined by continuity. 
• For (3 = 1 

S 1 ) in.. r,. \ — Df r^(l) 



Pn,2N( Xl ' 



^ i^2N x j)]7,j=l 



where the skew-symmetric 2x2 matrix kernel for x ^ y is given by 



a 



c. 



^(x-y) (D(x)D(y)) (i) 



l 

x—y 



D(x) 



c: 



(i) 

2JV+2 



a 



or 



i 

x-y 

-y) 



\D(0/ A W 

2N 



C=y _ 

(i) 

2JV+2. 



and for x = y the kernel is defined by continuity. 
• For (3 = A 

where the skew-symmetric 2x2 matrix kernel for x ^ y is given by 



K ( *\x,y) 



C ^( x -y)(D\x)D\y)). w 



2C 



N 



D*(x)\ 

flf)/ A ( 4 ) 



c 



2(^-2/) 

and /or x = y the kernel is defined by continuity. 
Proof. The Sokhotski-Plemelj formula implies that 



(4) 
JV+1 



2C 



(4) 



-v) 



D 2 (x)\ 

W)/ A (4) 

A 



C=IL 



C=x,v=y- 



PtiJn( Xi ' • • • ' x «) 



du\ ■ ■ ■ du r 

gn 



du\ ■ ■ ■ du r 



%)•••%) 

D 2 ( Ul )...D 2 (u n ) 
D*(vi) — D*(v n ) 



A 



(0) 



Vl=Xl, 



A 



(0) 



v 1 =x 1 ,...,v„ 



We use the first formula for /3 = 1,2 and the second one for (3 = 4, and then proceed similarly to 
Propositions EU and EHU □ 

Remark 4.3.2. If instead of asking for the measure [i to have nice density near x\, . . . ,x n , we 
require that fi is purely atomic near these points, then the formulas of Theorem 14.3.11 will continue 
to hold if we understand the symbol [/(C)]cc as the residue of the meromorphic function f(() at 
the point x. The proof is very similar. 



The one- and two-point averages 

1 



(D(x)) , 



(D(x)D(y)) , 



D(x) 



1 



M*)/' N V ~' \D(y)/' \D(x)D(y) 

for (3 = 1,2, and similar expressions with D 2 instead of D for (3 = 4, can all be evaluated in terms 
of (skew)-orthogonal polynomials associated with the measure /U. 

For (3 = 2 one needs the usual orthogonal polynomials associated with [i. The computation has 
been done by Strahov and Fyodorov [HH], see also pj. 

For (3 = 1 and (3 = 4 the skew-orthogonal polynomials associated to the skew-symmetric inner 
products 

(f,9)l = o [[ S S^( X ~ y)f( x )9(y)p(dx)fi(dy), (f,g) 4 = ■= [ (f(x)g'(x) - f'{x)g{x))fi(dx) 

1 J JR' 2 1 JR 

need to be constructed. The expression of the one- and two-point functions in terms of these 
polynomials are exactly the same as in (3 = 2 case modulo the change of the inner product. Since 
the proofs are essentially the same, we omit them here. 

If one substitutes the expressions in terms of the (skew)-orthogonal polynomials into the formulas 
for the correlation kernels in Theorem 14.3.11 one recovers the classical formulas of Random Matrix 
Theory, see e.g. 
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5. ASYMPTOTICS 

This Section presents the asymptotic analysis of kernels and correlation functions of characteris- 
tic polynomials. The three cases are considered: the case of the Gaussian Unitary ensemble (GUE), 
the case of the Gaussian Orthogonal Ensemble (GOE), and the case of the Gaussian Symplectic 
Ensemble (GSE). 

5.1. Statement of Asymptotic Results. 

5.1.1. The Gaussian Unitary Ensemble. In the notation of Section^ we consider the case of (3 = 2, 
fi(dx) = e~ x dx. The computations of correlation functions for GUE involve monic orthogonal 
polynomials defined by the following inner product: 



+oo 



</,<?>; 



dxe x f(x)g(x), 



and the relevant family of (monic) orthogonal polynomials {7r n }, n = 0, 1, 2, . . ., has the property 

Remark 5.1.1. The inner products as well as the averages over ensembles are denoted by the 
same symbol () throughout the paper. These notations should not lead to a confusion. 

In Section 0] we have found that averages of products and ratios of characteristic polynomials 
can be expressed in terms of kernels. In the case of GUE these kernels are 



1 



h 



(2) 



N 



W; 



C-V \D(v)/ A P)' 



w, 



(2) 



NJII 



(c, f; ) = e i (c-,)< i5 ( ? ^ M > a(a ,c, f; €c\ R 



The following statement was proved in Ref. |53) . 

Theorem 5.1.2. (Scaling limits of kernels at the origin of the spectrum). The following limit 
relations hold: 



hm -i=W#i( " 



w^oo V2iV 



= ^GUe(C,v), Km-Lw^ 



- BUI 



li m -^=wS ] - \rr(-F= . ~F=) = e(C ,V) 
n^oo^/2N N ' IH W2N\/2N gue ^> " 



where 



®gue(C,v) 
$gue(C,v) 



IT 



sin(C — t/) 



(e C, rje C. 



Iran > 0, 



f-^o - C e c, r, s 

Imi] < 0. 



§Gue(C,v) = 2™< 



AC-v) 

C-r/ 
e -i{C-v) 



Im C > 0, Im rj < 

Im C < 0, Im 77 > 
otherwise. 



C,v£ 



c 



v 



) = §gue(C,v), 
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Theorem 5.1.3. For any integers N > 1 and S > 1 — N, and finite sets of mutually distinct 
complex numbers 

a~ = {ajf ,. . a + = {a{ , . . ..a^}, 
f3- = {fc,...,[}- 2 }, /?+ = {/?+ ...,/?+}, 
with \a~\ — \a + \ = \/3~\ — \(3 + \ = S and 

a -C\a + = 0, a + nR = p + nM = 0, /Tn/3 + = 0, 



one has 



lim < 

AT— >oo 



T(2)(iV) \ fel ^/ajNrrfa Df 




where 



1 > 



T(2)(iV) M2) 



(2N)T, 



&gue{& , a+ ) is a matrix with rows parameterized by elements of a and (3 + , columns 

parameterized by (3~,a + , and with matrix elements: 

• 8 GUE (a-,(3j) =S I GUE (a~,f3-) 



• S GUE (ft,[37)=&J UE (p-,ft) 
. § G ^(/^a+)=§^(/?+ a+) 

The constant C N can be explicitly computed, see Mehta [H| : 

( N N(N-l) JL 



3=1 

5.1.2. The Gaussian Orthogonal Ensemble. Let us now consider the case (3 = 1, ^i(dx) = e 2 dx. 
The computations of correlation functions for GOE involve skew orthogonal polynomials. They 
are defined by the following skew symmetric inner product: 



-00 



(f,g)i = 2 I ° Ln ■ I s " u( w ~ ''•)/(•'').'/(//)• 

Specifically, the relevant family of (monic) skew orthogonal polynomials {p n }, n = 0,1,2,..., is 
given by 

(P2k,P2l+l)i = ~ {P2l+l,P2k)l = h^6 kh 

and by the condition that all other skew symmetric products between these polynomials are zeros. 

The averages of products and ratios of characteristic polynomials can be expressed in terms of 
kernels. In the case of GOE these kernels are 

• < 1} /(C V) = 4t (C - rj) (d(C)d(v)) A d) C G C, n G C. 

N 

• wjtUcr,) = -J- ; Cg c,,€ 

C-r? \<?7)/ A W 
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The kernels W^\((, 77), W^\j((, 77), and Wjp ni (£, rj) admit representations as Christoffel-Darboux 
type sums. Namely, 



(5.1.1) 
(5.1.2) 



w$AC,v) = £ 



N 



P2i+l(()P2i(v) - P2i(C)P2i+l(v) 



i=0 



(1) 



1 N ~ 1 



P2i+i(()h2i(v) - P2i(C)h 2 i+i(v) 



i=0 



h 



(1) 



(5.1-3) < ) m(C,r ? ) = E 



JV-2 



h2i+i(0h 2 i(v) ~ h 2 i(C)h2i+i(v) 



i=o 



(i) 



Here i?^ = (( — x) 1 , and the functions /ifc are the skew symmetric products between skew orthog- 
onal polynomials and R^: 

(5.1.4) h k (ri) = { Pk ,R n ) 1 , V € C\R. 

Theorem 5.1.4. (Scaling limits of kernels at the origin of the spectrum). The following limit 
relations hold: 



Km <>„(— — 



,ere 



S# OJ5 (C,»7) = <! 



2JV V2iV 
sin((" — ??) 



C-v 



C-r? 
Imr/ > 0, 

Imry < 0. 



C-r/ 
' +?° e^-^dt 



I , 

1 1 



Cg C, r,e C. 

Im( > 0, Imr/ < 



+00 p -i(C-v)tf]f (■ „ c 

- / ImC < 0,lmr? > ' ^ G 

1 * 



otherwise. 
Theorem 5.1.5. For any integers N > 1 and S > 1 



complex numbers 
such that 
one has 



N, and finite sets of mutually distinct 
a = {a 1 ,...,a k }, (3 = {ft, . . . , /3 m }, k - m = 2S 



pnw = 



lim 



1 / ntr^) \ 1 n(a;/3) 

T«(iv)\nr=i^)/ A(1 )j 



SGOs(a,/3|a,/3) 



where 



TW(N) c {1) 



CiVr km fc(fc+l) m(m-l) 

- A (27V) "2 4 4 

2Af+25 



SgoM ) /3|ck, /3) is a skew symmetric (k+m) x (k+m) matrix with rows and columns parameterized 
by elements of a and ft and with matrix elements: 
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if) 



§>GOE(.oci, otj) = § GOE (ai,atj) 

§GOE («» , Pj ) = SgO E( a i'Pj) 
S G OE(P l ,P j )=$ I G OE(A^j) 



The constants cijl can be explicitly computed, see Mehta 



(2iV)!C« = (20" J] 



j=0 L[ -2> 

5.1.3. T/ie Gaussian Symplectic Ensemble. Finally, let us take = 4, fi(dx) = e~ x dx. Given two 
functions f,g the skew symmetric inner product for the Gaussian Symplectic Ensemble (GSE) is 
defined by 



(/> 9) 4 



+oo 



(fee"* 2 (/(x)</(x) - . 



The family of monic skew orthogonal polynomials associated with GSE is introduced by 

{p 2 k,P 2 l+l) 4 = - (P2Z+l)P2fe)4 = ^fc^fcZ- 

All other brackets are zeros. In the case of GSE the kernels that determine averages of characteristic 
polynomials are: 



1 



h N 

l / D 2 (C)\ 

1 x ,(eC,r)e c\R. 



C-r? \L> 2 (r ? )/ A (4) 



<%(C^) = ^-i(C-^)( 



1 \ 



D\C)D\ri)J 



A 



(4) 



, C G C\M, 77 € 



Similar to the case of GOE the following expressions were obtained: 

N 

(5.1.5) 



I I'M i=0 



/I 



(4) 



(5.1.6) 



W 



P2i+l(()h 2 i(v) ~ P2i(()fl2i+l(v) 



i=0 



/J 



(4) 



(5.1.7) 0-l) = h ff-i 



,(4) 



■JV-2 

E 

i=0 



h2i+i(0h 2 i{v) ~ h 2 i{C,)h 2 i + i{ri) 



h 



(4) 



+ (R(, R 



Theorem 5.1.6. (Scaling limits of kernels at the origin of the spectrum). The following limit 
relations hold 



JV-oo 7V ' J W2N ^/2iV 



1 w(4) / C ?/ 



lim — WS } TTT (—=, 



,ere 



„r v 1 r , sin(C — n)x _ _ 

S^(C,»7) = -/^— — > CG C, 1)6 c. 



7T 



o 



r,0 



-III 

GSE 



(C,r?) = 2ni < 



d_ 

d 

dC 
o 
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3 -i(C-v) 1 

Im C < 0, Im rj > 

(~V J 

e i(C-??) 1 

Im C > 0, Im rj < 

C - V J 

otherwise. 



Theorem 5.1.7. For ant/ integers N > 1 and S" > 1 — iV, and finite sets of mutually distinct 
complex numbers 

{ai,..., a fc }, /3 = • • • ,/3 m }, k-m = 2S 



a 



such that 



one has 



/?n: 







lim < 



tW(jv) \n™id( 




n(«;/3) 

V(a)V(f3) 



Pf 



§G5£;(a,/3|a,/3) 



where 



C{9 km fc(fc-l) m(m+l) 



T( 4 )(iV) /-W 

§G5£;(a, /3|a, /5) is a skew symmetric (k+m) x (fc + m) matrix with rows and columns parameterized 
by elements of a and 0, and with matrix elements: 

• S GS E(ai,Uj) =S I GSE (a i ,a j ) 

• $GSE(Pi,P j )=S I G I £ E (l3 i ,l3 j ) 



The constants cff can be explicitly computed, see Mehta [ 

_ (2^)^+1/50 



iV!C 



2iV(JV+l/2) 



7V-1 
3=0 



The following Sections include the proofs of Theorems 15,1.41 and 15.1,61 

Remark 5.1.8. As we have seen in Sections 13 . 9114, 3l the correlation functions of GUE, GOE, GSE 
can be expressed through 2-point averages of characteristic polynomials, see Theorem 14.3.11 If 
we insert the asymptotic expressions for Sque, ^goe, ^gse to the formulas of Theorem 14 . 3 . 1 1 we 
obtain the following asymptotic formulas for the correlation kernels: 



lim 



1 



-.K 



(2), 
N * 



x y , 
/2N' V2N' 



sin(a; — y) 
ir(x - y) 



lim V2JV 

AT->-oo \ 




V 



d sin(z-y) 
dx x—y 

s\n(x—y) 
x-y 



s\n(x—y) 
x-y 

+oo 
r dt s'm(x—y)t 

J t 
1 
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The expression in the righthand sides are well known 3 , see e. g. Forrester |27| . Chapter 5. To 
derive these expressions we observe that 

+00 

l*>GOE&,V)\ <=Xtri=v - ~ J 1 

1 

PGSE{Qr,)\ c=x>r)=y - x _ y 
5.2. Summation formulae for kernels. 

5.2.1. Summation formulae for kernels in the GOE case. 

Proposition 5.2.1. The Christoffel-Darboux type sum constructed from the skew orthogonal poly- 
nomials {p n }, n = 0,1,2,..., is representable in terms of the orthonormal functions associated 
with the Hermite polynomials: 

^ P2i( X )P2i+l(y) ~ P2i(x)P2i+l(y) 

; 5.2. 1 : i=o i 



x 2 +y 2 Cat 



e 2 



CJV-1 



+ VN-i{y)VN{x) 



dx\ x — y 



where the functions i/j n (x) are orthonormal functions associated with the Hermite polynomials. 
Namely, 



- 2 £_ 

2 



(5.2.2) ipn(x) = 7r "^ e = Hn( f] e ) / Tp n (x)ip m (x)dx = 5 nm . 

Here H n (x) are the Hermite polynomials, 7r n are the monic orthogonal polynomials defined by the 
weight w(x) = e~ x , and the coefficients c n are their norms, 

(5.2.3) c n = tt 1 /^-"/ 2 ^!. 



This result is well-known and can be found in the literature on Random Matrix Theory (see, 
for example, Forrester, Honner and Nagao Adler, Forrester, Nagao, and van Moerbeke [Q, 
Widom jSH]). 

Proposition 15 . 2 . fl lead us to the representations for kernels which are suitable for the asymptotic 
analysis. 

Proposition 5.2.2. The kernels W^\((,r]), Wjpjj(£, rj), and wffl in (C, rj) are representable in 
terms of the orthonomal functions ip n ( x ) defined by equation \5.2.M\) . and by their Cauchy type 
transforms ^ n (C), 

(5.2.4) ».(0 = f ±£^HwW . Ce 



x — C 



Namely, 

(5.2.5) <(U) = -^e^ 

CN 



Tc{ J+^)W(C) 



3 The correlation functions in j3 — 1,4 cases are usually written as quaternion determinants while we use pfaffian 
representations. The relation between them is explained in Section [ 



52 



A. BORODIN AND E. STRAHOV 



CN_ e? 2 /2 



(5.2.6) 



CJV-1 



W 



(i) 



NJII 



(C,r?) = - 



1 C N - 



(5.2.7) 



2 CAT_2 
*iV-2(C) 



+ -^-i(C) 



^AT-l(C) 



C-r? 

dccdj/ 4^N{x)e~ y2 ^ 2 sgn(y — x) 



dxdy e~ y2 1 2 sgn(y — x)ipN-2{x) 



-1/2 



(y - »7)(a: - 
dxd?/ e~ y2 / 2 sgn(y — x)ipN-i(x) 

(y - v)(x - C) 

dxdy e~ y2 / 2 sgn(y — x)if)N-i(x) f f dxdy e~ y2 / 2 sgn(y — x)ipN-2(x) 



II 



(y-0 J J (y-v) 

Proof. Equation (|5.2,5|) is obtained immediately from equations 1)5.2.1(1 and (|5.1.1|) . Let us derive 
formula 1)5.2.6)1 . We start from the representation of the second kernel as the Christoffel-Darboux 
type sum, equation ()5.1.2)1 . In this representation we rewrite the functions hn ij]) , ^2i+i (v) explicitly 
as integral transformations (equation ()5.1.4|) ) of the monic skew orthogonal polynomials p2i and 
P2i+i- Then we apply formula ()5.2.1|) and obtain 



(5.2.8) 



C-r] 2 c N -i 
dxdyeT y2 l 2 sgn(y - x) d ( %1) N {x)^n-i{0 - ip N -i(x)ip N (() 



v-y 



dx 



x — ( 



dxdyi^N(x)e y2 1 2 sgn(y — x) 



v-y 

The first term in the brackets is simplified if we perform the integration by parts. After that we 
use 



1 



1 



1 



1 



{y - Q{y -v) C-v \y-C y-v, 

and decompose the resulting integral into two terms. In each term the integrals are the Cauchy 
type transforms ^f n of functions ip n . Now we apply the well-known formula 

fiv(C)Vw-i(C) - *jv-l(C)^iv(C) = — e~ cV2 - 

c N 

This formula follows, for example, from the fact that the determinant of the solution of the 
Riemann-Hilbert problem for orthogonal polynomials is identically 1, see Deift page 44. Then 
one term (the first term in equation ()5,2.8)0 is cancelled with and equation (|5.2.6|) is obtained. 
Equation (|5.2.7|) is derived from equations (|5. 1 .3|) and ()5.2.1)) by the same procedure. □ 



It is convenient to introduce two functions, 

In{v) ~- 



J J dxdy e y I 2 sgn(y — x)ipN-i{x) 



and 



E N (C,V) 



y-v 

dxdy e~ y2 1 2 sgn(y — x)ipN-i(x) 

(y-v)(x-Q 



Then the kernels W^ H (C,v) an d ^Ninid v) can be rewritten as follows, 



(5.2.9) 



w£Uc,v) 



- e s / 

CN-1 



NJII 

V N (v)lpN-l(0 ~ ^N-i(v)^n(C) 



C-r] 



+ ^n-i(()In(v) 
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and 

(5.2.10) 



2 CN-2 



En^Cv^nMO ~ E n _ 2 (C,v)^n-i(C) - -In-i(()In-2(v) 



5.2.2. Summation formulae for kernels in the GSE case. Given two functions, f{x) and g(x), 
introduce new skew symmetric brackets (,) by 

(f,g) = l [ dx(f(x)g\x)-f'(x)g(x)) 



(5.2.11) 
It is easy to check that 

(5.2.12) (f,g) 4 = {f{x)e- x2 /\g{x)e- x ' 2 / 2 ). 

Orthonormal functions associated with the skew orthogonal polynomials are defined by 



1 

e 2 p 2i (x), <j)2i+l{x) 



1 



h 



(-1) 



h) 



(4) 



e 2 p 2i+ i{x) 



Relation (|5.2.12() implies that 



dx (4>2i(x)(/) 2j+1 {x) - (j)2i + l{x)<t> 2 j{x)) = 5ij, 



and all other possible brackets are zeros. The following formula is well known, see Forrester, Nagao 
and Honner[2*nj. Widom 58,, and Adler, Forrester, Nagao, and van Moerbeke [Q 

N-i r 



(5.2.13) 



i=0 



4>2i{x)(p2i + i(y) - 4>2i+i(x)4>2i(y) 



C2N 
C2N-1 



x-y 



+ ^2N{y) I 1p2N-l(t)dt 



where the functions ip n {x) are orthonormal functions associated with the Hermite polynomials (see 
equation (|5.2.2|0 . 

Proposition 5.2.3. The summation formula for the first kernel, W^j(C,r]), is 
(5.2.14) 



C2N+2 



e 2 



C2JV+1 

c c 

1p2N+2(l])' t p2N+l(x) ~ 4>2N+l{ri)4>2N+2{x 



dx 



Tj — X 



+ / ^2N{x)dx I 1p2N-l(t)dt 



Proof. Application of basic summation formula (|5.2.1MJ) together with expression (|5,1.5j) for the 
kernel. □ 

Proposition 5.2.4. It is possible to rewrite the functions h2i(f]), /i2i+i( ? ?) as follows 



(5.2.15) 



X — T) 

dx 
x — r\ 

Proof. Apply relation (|5.2.12j) and integrate by parts. 



J 2i\ 



(5.2.16) 



h 2i+ M = ^ 



" a;2 /2j 1 / 



<t>2i+l( x ) 



□ 
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,(4) 



Proposition 5.2.5. The second kernel, W*, jj{CiV)> can be rewritten as follows 
(5.2.17) 



C 2 /2 C 2N 
C2N-1 



^2N(y)lp2N-l(C) ~ ^2jV-l(^)^2Af(C) 

C- v 



+ ^2N(V) / 4>2N-l{t)dt 



where ^> n was introduced by equation \5.2.1$ 

Proof. We insert the expressions for the functions h2i(n), /i2i+i(?7) obtained in Proposition I5.2.H 
into the formula for the kernel W^\j(C,rj), equation (j5.1.fi|) . After that we use 

1 1/1 1 



{y-(){y-v) C-r/\y-( y-vj' 

and decompose the first resulting integral into two parts. Then we apply the well-known formula 

(5.2.18) ^2N(0^2N-l(C) ~ *2AT-l(0^v(C) 

and derive equation (|5.2.17|) . □ 



C2N-1 g -C 2 /2 
C2N 



It is worth noting that the first term in the brackets of equation 1)5.2. 17jl is exactly the same as 
in the case of GOE, see (|5.2.9j) . 



Proposition 5.2.6. Let £ £ C\R, rj E C\R. The summation formula for the third kernel, 

7 W ( 
N,III\ 



(5.2.19) 



C2JV-3 



^2N-3(v) F 2N-2(riX) ~ ^2N-2(v) F 2N-3(ll,C) 
~ ^2N-2(V)^2N~3(C) , 

dx ipN{%) d e~ x2 l 2 



where we have introduced: 

(5.2.20) F N ( V ,() , , 

n — x dx Q — x 

Proof. Equation ()5.1.7|) shows that we need to compute the following kernel: 

A 2 h2i+i(()h 2 i(ri) - h 2 i(C)h2i+i(v) 
=o 



(5.2.21) 



^ / (C,r ? )=E 



h 



(4) 



We rewrite the functions h 2 i+i((), h 2 i((), h 2 i+i(r]), h 2 i(r)) explicitly using equation ()5.2.15|1 and 



d e-^l 2 
' dx C, — x 



h 2 i+i (C) = V hf ] I dx 4> 2l+ i (x) 



d e~ x2 /\ 
dx C — x 



dx 
x — r] 



dx 
x — r] 



,J 2i+l 



X . 



Insert the above expressions to equation (|5.2.21|) . and obtain 



N-2 



s i n ii (c^) = Y1 dxdy 



i=0 



e -yi/2 d e - X */2 
rj — y dx ( — x 



{^2i+l(y)^2i{x) - (/>2i(y)<hi+l{x)) 
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Now we apply basic summation formula (|5.2.13|) and find 
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C2N-2 
C2N-3 



dxdy 



e-y 2 ' 2 d e~ x2 / 2 



rj — y dx Q — x 



+ 



dxdy 



1p2N-2(x)lp2N-3(y) ~ ^2N-3(x)lp2N-2(y) 

x-y 

x 

i>2N-2{y) / 1p2N-3(t)dt 



e~y 2 l 2 d e~ x2 l 2 
rj — y dx C — x 



C2N-2 
C2N-3 



Rewrite h(C, v) 

repeated integral 

d e- 2 /2 



(5.2.22) h(C, V ) 
Apply identity 



dx 



dx C — x 
1 



dy 



-y 2 / 2 1p2N-2{x)lp2N-3(y) ~ tp2N-3{x)lp2N^2 (y) 



v-y 
l 



x-y 



l 



l 



iv-y)(x-y) rj-x\x-y r]-y / 
in the inner integral of equation (|5.2.22j) . Then we obtain, using equation (|5.2.18ft . 

e^y 2 ! 2 1p2N-2(x)l^2N-3(y) ~ 1p2N-3(x)lp2N-2(y) 



dy 



1 



1] — X 

1 

77 — X 



dy e 
dy e 



v-y x-y 

y 2/2 ^2N-2{x)ll>2N-3{y) ~ ^2N -3{x)^2N -2{y) 

x-y 



1 



rj — x 



-e-* 2 / 2 + 



1 



V/2 Tp2N-2{x)lp2N-3{y) ~ 1p2N -z{x)lp2N -2{y) 

v-y 

(4>2N-2(x)$>2N-3(v) ~ 1p2N-3(x)^2N-2(v)) 



Therefore, 



h(C,v) 



+ 



dx 



rj — x 



d e~ x2 l 2 1 



dx C — x v — x 
d e- x2 ' 2 1 
dx Q — x rj — x 



e~* 2 / 2 + 



{i'2N-2{x)^2N-3{v) ~ 1p2N-3(x)^2N-2(v)) 



or 



h((,v) 



1 



1 



+ *2JV-3(?7)-^2JV-2(r?, C) - ^2N-2(V) F 2N-3(V, C) 

t V — x Q — x I 

We simplify /2(C) V) by partial integration, obtain an expression for 5^/ 7 (C, rj), insert this to formula 
(15X71) . and obtain (15.2. 1911 . 

5.3. Asymptotics of functions associated with Hermite polynomials. 



□ 



5.3.1. Asymptotics of ip]y(x). The large N asymptotics of the functions ipN(x) in different regions 
of the real line can be found in Szego [S3]) formula (8.22.14). However, the asymtotic expression 
near the edges of the spectrum is valid only in infinitesimally small neighborhoods (as N —* 00) of 
the edges. In what follows we need uniform asymptotics of functions tpjsr(x) on the entire real line. 
Theorem 2.2 in Deift, Kriecherbauer, McLaughlin and Zhou 23 provides the needed asymptotics 
of orthogonal polynomials with respect to the weights w(x)dx = e~®( x ^dx on the real line, where 
Q(.x) = X^fc=o^ xfc ' ^2m > (see also Deift, Kriecherbauer, McLaughlin and Zhou |H]). The 
expressions in Theorem 2.2 are valid in the entire complex plane. We adopt this result to the 
particular case of Q{x) = x 2 . 

The parity of functions tPn(x) defined by equation (|5,2.2|) corresponds to the parity of N, i. e. if 
N is even, ip^{x) is an even function, and if N is odd, iPn(x) is an odd function. Thus it is enough 
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to consider the asymptotics of functions iPn(x) on R + . We decompose M + into three regions: B$, Cs 



B& 



Cs 



As 



1-5 1 1 + 6 



Fig. 4. Different asymptotic regions for ipN(y/2Nx) 



and As, see Fig. 4. 

Theorem 5.3.1. (Plancherel-Rotach asymptotics of i[)n(\/2Nx) on the real line.) There exists a 
So such that for all < 5 < 5o the following asymptotic formulae hold: 
(i) For x £ As 



(5.3.1) 

(ii) For x G Bs 

(5.3.2) 

ip N (V2Nx) 

( Hi ) For x G Cs 



ij N (V2Nx) 



23/4^1/2^1/4 

ttN 



x-l 



1/4 



+ 



X + 1 



1/4- 



x exp 



x + 1 j \x — 1 

\/x 2 — 1 — ln(x + \J x 1 — 1) 



l + ^(-) 



: 7r i/ 2 2 ^i/4 (1 _ ^2)1/4 { cos (^V 1 - x 2 + (iV + 1/2) arcsin x - ^) 
x ^1 + C(-^r)) + sin (Vx Vl -x 2 + (N- 1/2) arcsin x - 0{^) 



MVwx) =j^\^0^Un{ X )) 1/a ^Un{x)) (i + o(^)) 



(5.3.3) 



(x-l 



,1/4 



(x + 1) i/4 (/iv(-))- 1/4 Ai'(/iv(x))^ + 0(-; 



where the function /at(x) is given 6y 



(5.3.4) 



37V \ 



2/3 



X 



\/x 2 — 1 — ln(x + vx 2 — 1) 



and Ai denotes the Airy function determined as the solution of 
(5.3.5) Ai"(x) = xAi(x), 

satisfying 



2/3 



(5.3.6) 



lim Ai(x)\/ivrx 1 / 4 e 



2 3/2 



1. 
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5V 



All the error terms are uniform for 5 £ compact subsets of (0, So] and for x G Xg, where X = 
A,B,C. 

Remark 5.3.2. When \x — 1| < e, e is a small parameter, /at(x) = (x — 1)(</>jv(x)) 2//3 , where 4>n{x) 
is a strictly positive function at x = 1. Therefore the function i^nW^Nx) is finite at x = 1. 

Corollary 5.3.3. sup\^n(x)\ < const -iV~i2. 

Proof. This is evident from the asymptotic expressions above for i/jn(x), and from the asymptotic 
expansion of the Airy function and its derivative. □ 

M 



Proposition 5.3.4. Assume that x £ [-M, M]; M — ► oo as N — > oo, 



as N — > oo . T/ien 



(5.3.7) 



ip N (x) 



2 l/4 



-I^TT-f)^(S) 



^1/2^1/4 

The error term is uniform for x £ [—M,M]. 

Proof. This formula follows directly from equation (j5,3.2j) , where we replace x by 



>2N 



in powers of 



'2N 



the fraction 



1 x z 

1 IN 



1/4 



, and the expression inside the cosine. 



and expand 



□ 



5.3.2. Asymptotics o/^jv(C)- Recall that the Cauchy type transform ^n(C) of the function iPn{x) 
was defined by equation (|5.2.4|) . 

Proposition 5.3.5. Assume that ( £ C\R, M — > oo as N -> oo, In N < M < v^/V as iV -»• oo. 
T/ien i/ie following asymptotic expression for \&jv(C) holds: 



(5.3.8) 



*iv(C) 



2 l/4 



7T 



l/2jyl/4 



e 2:2 / 2 



cos 



2 



M 3 In M 

iV 3 / 4 



TTte error term is uniform for all £ taken from a compact subset of C\M. 



Proof. We decompose the integral in the definition of the function \I/jv(C) (equation (|5.2.4j0 into 
three parts: 



M 



M 



(5.3.9) 



dx e x / 2 ip n (x) f dx e x / 2 ip n (x) f dx e x / 2 ip n (x) f dx e x / 2 ip n (x) 



x-C 



-M 



X — ( 



+ 



M 



x-C 



+ 



By Corollary 15.3.31 ibiv (x) is uniformly bounded by const -N 1 / 12 . This enables us to estimate the 
second and the third integrals in the expression above. Namely, 

+00 2 

dx e~ x > 2 ip n (x) 



I 

M 



x-C 



< const -AT- 1 / 12 



+oo 



dxe x / 2 : 



M 



const -iV-i/i2 e -MV2_ 



This shows that the second and third integrals in ()5.3.9|) are small. To obtain formula Q5.3.8|) we 
insert the asymptotic formula for V'aK 3 -) (equation Q5.3.2|) ) into the first integral in (|5.3.9j) . and 
note that 



M 



(5.3.10) 



-M 



dx e x2 l 2 



M 



< 



-M 



dx 

W^c\ 



O(lnM) 



□ 
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Proposition 5.3.6. Under the same assumptions as in Proposition ^ 3. d\ the asymptotic formula 
for the Cauchy type transform vE'jv(C) of 'ipN(x) is 

C v 2 1 / 4 tt 1 / 2 \ie^-^\ ImC>0 fMHnM 
(5.3.11) * N (-=) = —. -. \ . ■ (t _ E N, ». n +£> 



/ 2N J N 1 / 4 |_i e -i(C-V), ImC<0 V ^ 3/4 

The error term is uniform for all £ taken from a compact subset of C\M. 

x 

I 'roof. We replace x by _ ^jy m ^ e ^ e S Is ^ i n equation (|5.3.8[) . The obtained integral is finite as 

N — > oo and has the following limiting value 

+oo 



/ 



dx 

cos 



ttN 
X ~ — 



x-C 

— oo 

We compute this integral by the standard residue calculations and get formula ([5.3.11|) , □ 
5.3.3. Asymptotics of Taj ( — = ) The aim of this subsection is to prove the following Proposition: 

Proposition 5.3.7. Let r] be taken from a compact subset K o/C\R. Then the following estimate 
is valid 

sup|/tv(^2=)| < const -N~ 1/2 
K V2N 

Proof. The function In {of) can be rewritten as 

r-+oo 

In(v) = / dxifijy (a;) cf)^ (x ) , 



where 



/+oo 2 
dye~ y / 2 sgn(y — x) 
-oo 5?- 



2/ 



M 

Let M — ► oo as N — > oo, — = as iV -> oo. Then the integral in the expression for In(v) can 
be restricted to the domain [—M,M], with an exponentially small error: 

cM 



(5.3.12) In(v)= dxV>iv(x)^(x) + 0(e- JU ). 

In order to show this we note that the integral can be rewritten as follows 

/+oo /"+00 
dye~ y /2 / sgn(y - x)ip N (x)dx. 

-oo fl-V J-oo 

We define 

r+oc 

etp N (y) = / sgn(y - x)ip N (x)dx. 



Decompose the integral In(v) into three parts: 
(5.3.13) 



M i p — M 1 

1 , , x / , _„2/o 1 



In(v) = / */e y /2 eipN(y) + / dye y /2 e^jv(y) 

»7-2/ 7-oo 

/H-oc 

+ / dye - * / 2 e^ N {y). 

Jm v-v 

The global estimate for e^>jv(y) is known (see, Nagao and Wadati [15] . equation (A.l)): 

sup |c0iv(2/) | < const • iV -1 / 4 
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r,9 



We use this estimate to evaluate the second and third integrals in equation Q5.3.13|) . For example, 
for the third integral we obtain : 

+C ° dye-y 2 ' 2 ^-e^ N (y) < ^ f°° e^ydy = 0{e~ M ). 
m rj-y iW 4 J M 

The same estimate is valid for the second integral in equation (J5.3.13|) . By equation Q5.3.12JI the 
problem of evaluation of Jjyf =) is reduced to the evaluation of the following integral: 



'2N' 



lN,M 



(77) = / dxip N (x)</> v (x), 

J-M 



We change the variables of the integration, x = t/y/2N, y = s/y2N, and find 



MV2N 



lNM 7m ] 



1 



*2N 



dti/> N (- 



f 2N' 



+00 

I dse- s2 / m sgn(s-t) — 

J v- 



-MV2N 

The last integral can be rewritten as 

+00 



MV2N 



lNM 7m ] 



1 



<2N 



dse- s2 ' AN - 



Tj — S 



dtipN(-f=)sgn(s-t). 



Now we can insert the asymptotic expression for ip 



MV2N 



f t 2 1 / 4 f ( 

J dt ^ N< "^/W^ Sgn ^ ~ t " > = 7^/2^1/4 J dtsga(s-t)i COS 



N 



MV2N 



-MV2N 
t 

(2iV)V2 



into the second integral and obtain 



irN' 



+ 



Af 3 
\/2N 



-MV2N 

This shows that the integral 



-MV2N 



MV2N 



dtipN{-^=)sga(s-t) 



-MV2N 

+00 



const -M ~"~ ° 2 1 
is bounded by —j- — . The integral J dse~ s l iN is bounded by const TogiV uniformly 

for all r\ £ K, where K is a compact subset of C\R. Choosing M = M{N) appropriately (such 



that 



— ^}}° gN < 1 as N — ► 00) we prove the Proposition. 

iI>n(x) dx 



□ 



5.3.4. Estimate of f 

M ^ 

Proposition 5.3.8. Assume that M <C \^N. Then the following estimate is valid: 

const const M 



f ?Pn{x) dx 


< Max 


J x 





M 2 jyl/12 ' N 3/4 



where the constants do not depend on N,M. 



Proof. We use representation of i^n(x) in terms of the Hermite polynomial H^(x) (see equation 
Q5.2.2J1 ). the relation H' N (x) = 2NHn-i(x), and take the integral j by parts. As a 



M 



X 



(,0 
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result we obtain 

+00 

i/jn{x) dx 



I 

M 



CN+1 



X 



(N + l)c N 



M 2 



+00 



+00 



ip N+ i(x) dx 



M 



M 



Since |V>at+i(M)| < const -N 4 / 12 f Corollary I5.3.3|) the first term in the brackets is estimated by 
const ■N~ l / 12 M~ 2 . As for the second term we note that 



M 



ip N+ i(x)dx 



M 



< 



ipN + i(x)dx 



+ 



ijj N+ i{x)dx 



< const -N- l,A + const -iV~ 1/4 M < const -N~ 1/4 M, 
where we have used the estimate proved in Nagao an Wadati |46j . equation (A.l): 



+00 



ij)N{x)dx 







< const -N~ 1/4 . 



The third integral is evaluated as follows: 

+00 



if) N+ x{x) dx 



x- 



M 



< 



M 2 



\ip N+ i(x)\ dx < 



const V2N 
N l / l2 M 2 ' 



M 



Noting that the factor 



cjy+i 



Proposition. 

5.3.5. Asymptotics of E]\r(C,r]). 
Proposition 5.3.9. Take (, rj G C\R. Then 



v , , , — 7==, as N — > 00, we obtain the estimate in the statement of the 

(N+l)c N V2JV' 



□ 



En 



(5.3.14) 



C V 



'2N V2N 



2 1 / 4 f dy e ~ y2/AN f dx cos ( x o~) sgn ( y ~ x 



Tri/ajVi/4 J Q- y 

' M 3 log 2 N losN 



r\ — x 



+ O Max 



as N — > 00. T/ie error term, O (Max 
compact subset of 



N 3/4 ' m 2 N 1 / 12 
M 3 log 2 N log N 



N 3/4 > MZtfL/Vi 



Remark 5.3.10. The error term O I Max 
in comparison with the main term, 



M 3 log 2 N logiV 

N 3/4 > M 2 N \/\1 



is uniform for (, rj taken from a 
in Proposition 15.3.91 is negligible 



vriV\ 



2 1 / 4 f dy e~ y2 / AN f dx cos ( x o~) sgn ( y ~ x 



7T 1 



1/2JVV4 7 ( - y 



r\ — x 



M 3 log 2 A 1 logiV 1 

< 7T. & , < 



if 

( 5 - 3 - 15 ) N 3/4 - N i/V M 2 N X I X2 ~ AV4 

Once M = M(N) is chosen such that A^ 1 / 12 (log iV) 1 / 2 < M < iV 1 / 6 (log iV)~ 2 / 3 the conditions 
above (equation (|5.3.15jl ) are satisfied. 
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Proof. It is convenient to rewrite En 



c 



v 



f 2N y/2N 



as 



(5.3.16) 



En 



c 



V 



'2N V2N 



-y 2 /4N 



9n(; 



u 



I) 



'2N V2N 



dy 



C-y 



where 



(5.3.17) 



9n(- 



y 



'2N V2N' 



sgn(y - x)i/> N ( 



'2N- 



dx 



X — 7] 



Equation Q5.3.14[) means that the function ipN-i(x) in expression (|5.3.16|) can be replaced by its 
asymptotics in the center of the spectrum if we are interested in the scaling limit of En(CiV)- I n 
order to prove that such replacement is allowed we proceed as follows. First we will show that 
the integration over y in equation (J5.3.16|) can be resticted to a domain on which the asymptotics 
of ipN-i{x) is given by Proposition 15.3.41 Second we will observe that the integration over x in 
equation (|5.3.17f) can be restricted to the same domain as well. 

Let K denote a compact subset of C\R. The function gN{y, —7=) nas the following bound 



which is uniform in y € R and in r/ £ K, 



f 2N 



(5.3.18) 



sup< 5Ar(y, , ) \ < const -N a , a — real positive. 
t/skI y2N 



In order to prove this inequality we rewrite the integral in equation H5.3.17|) as a sum of three 
integrals: 



/2N(l+e) 



9n(v, 



I2N' 



sgn(y — x)iPn(x) dx 



+00 



+ 



(5.3.19) 



-V2JV(l+e) 

-V2N(l+e) 



/2N 



sgn(y — x)iPn(x) dx 



/2AT(l+e) 



/2N 



sgn(y — x)i/jn(x) dx 



I2N 



where e > 0. The following estimate is evident from the asymptotic expressions for the function 
i/jn(x) (see Section 15.3.1(1 . 

\^ N (x)\<e~ coiastVN ^, 4i>l + e,e>0. 

V2N 



Using this estimate we obtain that the second and the third integral in equation (|5.3.19|) are 
exponentially small. As for the first integral in equation (|5.3.19f> it is clearly polynomially bounded 
since \x - ^|=| > ^^=, (Imr] / 0), and sup|^Ar(x)| < ^°°f 2 (see Corollary 15.3. 3j) . Thus inequality 

(|5.3.18|) is proved. 

Let M — ► 00 as N — > 00, but 



M 



as N — > 00. Then we claim that 



(,2 
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where the error term, 0(N a e M V 2 ) ; is uniform for £,77 taken from a compact subset of C\R. 

C V , 

as 



Indeed, let us take M as described above and rewrite En 



'2N V2N 



En 



c 



V 



'2N V2N 



M 



-y 2 /2 



9n(v, 



/2N- 



dy 



-M 



+ 



/2N 



-M 



I2N 



+ 



M 



I2N 



'- y2/2 9N(y, T ^)dy 



and estimate the second and the third integrals. Namely, 



Te- y2/2 9N(y,^)dy 



M 



'2N 



< 



+00 



~ y2/2 9N(y, ~^)y d y 



M 



+00 

< const -N a J e~ y2/2 y dy = const -N a ■ e" A/2/2 . 



where equation (j5.3.18j) was used. This shows that the integral over y in the expression for 

En(-^^, can be restricted to the domain [-M, M] with a small error term. 

Now assume that y is taken from the interval [— M, M], and M = M(N) is chosen such that 

M 

M — > 00 as N — > 00, but —^= as JV -> 00. Then 



N 



(5.3.20) 



sgn(y — x)^n(x) dx 



At 



/2N 



sgn(y — x)iPn(x) dx 



-M 



/2N 



+ 0{ Max[MA^" 3/4 , M~ 2 N~ 1/12 ] 



as N 00. The error term, ojMaxfMN-^M^iV- 1 / 12 ] j, is uniform on y £ R, r/ E C\R, 

where 77 is taken from a compact subset of C\M . In order to prove (|5,3.2Uj) we decompose the inte- 
gral in the left-hand side of equation (|5,3.2Uj) into three integrals on intervals [— M, M], (—00, — M), 
and (M, +00). The integrals on intervals (—00, —M) and (M, +00) are evaluated using Proposition 

EMI 

To summarize we have proved that 



E N {- 



c 



V 



'2N V2N' 




sgn(y — x)ipN(x)dx 



+ 



-M V2N ~ V-.W 

-M 2 /2\ 



I2N 



Max[MA^- 3/4 , M- 2 N- 1/12 } 



+ 0{N a e 



■)■ 
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The error term 0(N a e M V 2 ) i s exponentially small. Integrating the error term O Max[ Jf /4 , M2 ^ 1/12 
we obtain: 

M 



V 



O 



M V2N 
AIV2N 



y 



MaxfMiV" 3 / 4 , M^N- 1 / 12 ] 



-y 2 i iN d y 
C-y 



o 



Max[MiY" 3/4 , M- 2 N- 1/12 \ 



< 



< o 



Max[M log NN- 3 ^, M~ 2 log NN- 1 ^ 12 ] 



Another error term is due to the replacement of i/jn(x) by its asymptotics taken from Proposition 
15.1141 This error term is evaluated as follows: 



MV2N 



-y 2 / 4N d y 
C-y 



Mv2N 



dx 
x — rj 



O 



M 3 



< O 



M 3 (log AT) 2 



AT3/4 



-MV2N 



-MV2N 

This completes the proof. 
5.3.6. Asymptotics of Fn(C,v)- 

Proposition 5.3.11. Let £,77 be chosen from a compact subset o/C\l 

F N (-^=,-±=) = 2N - 
(5.3.21) 2 ^ 2 ^ 



□ 



M 



n l/2 N l/4 



, TTiV. 

at cos(i —J 



— A/ 

+ O [Max 



Then 

1 d e~ t2 l iN 
rj-tdt C~t 

1 



where M = M(N) is chosen such that M(N) -^00 as N — > oo ; 



M 4 

Af(JV) 



as iV -> 00. 



Proof. Recall that i*jv(C> ??) was defined by an integral in equation ([5.2.20)1 . Change the variable of 
integration in this integral, and decompose it into three integrals. The first one runs from —00 to 
—M, the second one runs from — M to M, and the third one runs from M to +00. In the second 
integral replace the function ifj^^x) by its asymptotics in the bulk of the spectrum (see Proposition 

M 4 

15.3.4)) . The integration of the error term gives 0( — 777). The first and the third integrals are of 

1 1 d e- t2 / m 
order 0{ - j9 Ar1 /19 ). It can be seen rewriting the expression 1~~ — ~ explicitly, and using 



'M 2 Ny 12 '' " " ° r "~ rj-tdt (-t 

the fact that sup 

|V>jv(*)| < const -N~ 1 / 12 (Corollary 



□ 



Remark 5.3.12. The error term, O ( Max 



M 4 



1 



jV 3 / 4 ' jyi/ 12 M 2 

the first term in the brackets of equation ([5.3.21)1 . if the following conditions are satisfied: 



is negligible in comparison with 



M 4 



1 



1 



N 3/4 ^ N l/4 ' jyl/12 M 2 ^ N l/4 ' 

The first condition implies that M <C iV 1 / 8 , and the second equation implies that A 1 / 12 -C M. 
Therefore if M(N) is chosen such that A 1 / 12 < M < iV 1/s the error term is negligible in compar- 
ison with the main term. 



5.4. Asymptotics of the kernels: GOE case. 
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5.4.1. Asymptotics of Wff\(—^ 
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.<2N V2N 

Theorem 5.4.1. (Scaling limit of the first kernel) For £ G C, r] G C, 



lim , 

2V-*oc 2iV 



1 1 d 


sin(C — 77) 




. C - v 



'2N V2N' 

Proof. Use formula (|5.2.5j) . the following relation which follows immediately from the properties of 
the Hermite polynomials, 

together with the asymptotic formula for the functions t[) n (x), Proposition 15.3.51 



□ 



c 



v 



5.4.2. Asymptotics ofW^\j( 



/2N V2N 

Proposition 5.4.2. Take C,rj G C\R. Then 



lim < V2N 



C-r? 



I C- v 



J(v-0 



Imi] > 0, 



e -i(v-0, Im?7<0 



Proof. Apply Propositions 15.3^11 and 15.3.61 

Theorem 5.4.3. (Scaling limit of the second kernel) Let ( £ C, r/ £ C\R. Then 



□ 



lim —^Wi 1] < — 



V 



>2N V2N' 



e *(ij-C) 

C — v 

e -i(v-0 



Imrj > 0, 
Imrj < 0. 



Proof. This follows immediately from equation Q5.2.9|) , Propositions 15.4.21 15.3.71 and from the fact 

□ 



V 



^ + ° N N AT 

that ~ 4/ — as JS — > 00. 

cjv-i V 2 

5.4.3. Asymptotics of j j' - 



2N V2N 

Theorem 5.4.4. (Scaling limit of the third kernel) Take £ G 
(5.4.1) 

f +?° e^-^dt 



(i) 



c 



lim U' Y , , , , , , 



2iri < 



Im C > 0, Im 77 < 



. Then 



+00 e -i(C-v)tdt C G 

/ ImC < 0,lmr? > ' s 

1 * 

otherwise. 



77 G 



Proof. Recall that the kernel ^^///(C) 7 ?) is given by equation (|5.2,lUj) . In order to determine the 
large N asymptotics of Wn\ h (C, 77) consider the expression: 

^Ar(C,r?) = Siv-i(C^)*iV- 2 (C) - ^7V~ 2 (C,r?)^iv-i(C). 
The asymptotics for £?at(Cj r ?) (equation (|5.3.14j) ). and the asymptotics for ^at(C) (equation l|5.3.11|) ) 
give us the leading term of -Ajv(Cj v) : 



(5.4.2) 



2 f . dy e 



N f 



-y 2 /AN p i{C-x) 



I 



sgn(y - x) 



rj-y 
2 n dye-y 2 l m 



C-x 



I 



rj-y 



I 



-i(C-x) 



sgn(y - x) 



C — x 



ImC > 



Im C < 
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In what follows the integral representation for sgn(y — x) will play a role, 

e i{y-x)t dt 



sgn(y 



1 

in 



t 



ft 



If we insert this expression to equation (|5.4.2[) . and perform the computations formally (i. e. 



changing the order of integrations, and replacing e 



-y 2 /4N 



by 1) we obtain: 



(5.4.3) 



A N (- 



c 



v 



'2N V2N' 



Since the factor ~ CN 1 ~ as N — > 
large limit in comparison with An( 



~ 27ri < 



/2iV 



<2N 



1 * 







oo, and the term -Zjv-i 



ImC > 0, Im rj < 

ImC < OJmry > 
otherwise. 



(see Proposition 15.3.7]) . formula ([5.4. lj) follows. 



2N> V2N' 

However these computations involve the integrals which are not absolutely convergent, so formula 
(|5.4.3f) must be justified. In order to prove this formula rigorously we observe that the following 
holds: 



(5.4.4) 



sgn(y - x) 



C 



-i(y-C)t d t 



ImC > 



Indeed, the integral in the left-hand side of equation (|5. 4.4(1 is equal to 2 dx — - (since 



0, when Im<^ > 0). To check that J i 



i 'J 

-oo 



-< — "'— J » -< isequalto f ' 

observe that the derivatives with respect to y are equal to each other. Second, it is evident that 
the first integral tends to zero as y — * +oo. The fact that the second integral tends to zero when 
y — > +oo can be proved by partial integration. Thus it remains to compute the large N limit of 
the expression 



-i(y-Ot dt 



we first 



(5.4.5) 



dy e 



-y 2 /4N 



+00 



rj-y 



<v-()tdt 



ImC > 0. 



The integrand, as a function in (y,t), decays exponentially as max(|y|, \t\) — > oo Therefore we are 
allowed to change the order of integration, and expression (|5.4.5|) is equal to 



(5.4.6) 



e^dt f dy e-^-y 2 /^ 



v-y 



ImC > 0. 



The large N limit of the inner integral can be computed as follows. Take M = M(N) such that 
M(N) —> oo as N — > oo, but M 2 (N)/N — > as N — > oo. Decompose the inner integral into three 
integrals: the first one runs from — oo to —M, the second one runs from — M to M, and the third 
one runs from M to +oo. Then the first and the third integrals are of order O(-h)- As for the 
second integral, it asymptotically equals to (2ni)e~' irit , when Imr] < 0, and to zero, when Imr/ > 0. 
It can be observed considering the semicircle contour Cm U [—M,M], where Cm is a semicircle 
from the point M to the point — M passing through lower half of the complex plane. The Jordan 
Lemma implies that 



(5.4.7) 



dz- 



-z 2 /4N -izt 



< 



const e 



-M 2 /4N 



M -t 
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Since t > 1 we obtain that 
(5.4.8) 



dz- 



-z 2 /4N p -izt 



M 



□ 



'C M V ~ Z 

where the estimate is uniform for t E [l,+oo). Thus the inner integral in equation (|5.4.Hj) is 
determined by the residue at r] of the integrand (uniformly for t £ [l,+oo) with error of order 

O ( — J ). Replacing the inner integral by this residue we again arrive to formula (|5.4,3|) 
\M J 

5.5. Asymptotics of the kernels: GSE case. 

5.5.1. Asymptotic of the first kernel. 

Theorem 5.5.1. (Scaling limit of the first kernel) For complex £,r/ 



(5.5.1) 



c 



f] 



lim w£\( _. _, 



dt 



sin(C — n)t 



Proof. Propositions 15.2^31 gives exact formula for the kernel (equation (|5.2.14|) ). The second 

term in the brackets in equation (|5.2.14j) is the product of two integrals. The first integral is of order 
C>(./V _3 / 4 ), since ip2N( 2 ^j ) < const -iV -1 / 4 uniformly for u chosen from a compact real interval. 

(Note that the asymptotics of functions iJjn in a real neighborhood of zero can be extended to a 
complex neighborhood of zero , see Szego jHS]; Deift, Kriecherbauer, McLaughlin, Venakides and 

Zhou |23| . such that the bound ^ 2 jv( 2 J^y ) < const ■N~ l l 4: remains valid). The second integral is 

bounded by const -TV -1 / 4 (see Nagao and Wadati [3^1) equation (A.l)). Thus the second term in 



the brackets in equation ()5.2.14|) (where C is replaced by 



c 



, and rj is replaced by 



is of 



2VN 2VN 
order 0(N~ 1 ). As for the first term in the brackets of equation (|5.2.14jl it is of order 0(iV _1 / 2 ). 
Indeed, this term is an integral over a compact interval, so the functions tp2N+2, 1P2N+1 can be 
replaced by their asymptotic expressions ( Proposition I5.3.4j) in the bulk of the spectrum. Thus 
we observe that the second term in the brackets in equation (|5.2.14|) is negligible with respect to 
the first term, and formula (|5.5.1|) is obtained. □ 

5.5.2. Asymptotics of the second kernel. 

Theorem 5.5.2. (Scaling limit of the second kernel) Let £ £ C. Then 

Imrj > 0, 



( e i(V~0 



(5.5.2) 



1 



lim 

N->oc2VN 



W 



("1) 



c 



N,II\ 



2VN 2VN 



Kv-0 



C — 17 



Inn? < 0. 



Proof. The exact expression for the second kernel in the case of GSE is essentially the same as in 
the case of GOE. Thus the asymptotics is obtained by the same method. □ 

5.5.3. Asymptotics of the third kernel. 

Theorem 5.5.3. (Scaling limit of the third kernel) Let r] are chosen from a compact subset of 
C\R. Then 

' . d f f- iC - 

C v 



(5.5.3) 



lim 

N— >oo 



2VN 2VN 



2iri < 



d( \r] 
d 



C 



Im C < 0, Im rj > 



-ir\ _ 



Im( > 0,lmry < ' 
otherwise. 
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Proof. Consider equation (|5.2.19j) . Replace the functions ^2N-3, F2N-2, F2N-3 by their 

asymptotics given by equations (j5.3.11|) and (j5.3.21|) . The first two terms in the brackets of 
equation (|5.2,19|) are then represented by an integral which can be directly computed by residue 
calculations. As for the third term in the brackets of equation (j5.2.19j) it is asymptotically negligible 
in comparison with two first terms. □ 

Remark 5.5.4. The asymptotics of the kernels W^\(a, (3), wjy\j(a, (3), W^\ n (a, 0) is deter- 
mined in Ref. [HE]- Note that this asymptotics can be obtained exploiting the asymptotic expres- 
sions for ipN ( Proposition 15.3.4]) and for (Proposition 15.3.6]) . 

6. Appendix. Some useful facts from the Linear Algebra 

6.1. The formula for the Cauchy determinant. Let A = (a\, ... , a^), B = (pi, . . . , be two 
non-intersecting sets. Then the following formula is valid: 

(6.1.1) detf^-V' -H^WH) 



6.2. Determinant of the block matrix. If the matrix H has the following block structure, 

(6.2.1) H = ( I A Q 
with square matrices A and B, then 

(6.2.2) det H = (-)I A H B I det A ■ det B. 

6.3. Minors of the inverse matrix. Assume that two N x N matrices A and B are inverse to 
each other, 

(6.3.1) A = B~ 1 . 

Then an arbitrary minor of the matrix A corresponding to the rows a\,...,a r and columns 
f3\, . . . , (3 r can be expressed as follows: 

(6.3.2) det A(au...,ar\0i,...,Pr) = 



. det B (1, ... Ji, ... J r , ...,N\1, ... ,d u ... ,d r , ...,N) 

= ( — r =1 - 

v ' detB 

Here k denotes the k^ row (or column) which is removed from the matrix. 

6.4. The expansion of the minors of the matrix I+A. Let a = (di, . . . , d r ), = (/3i, . . . , /? r ) 
denote removed columns or rows. Assume that a D $ is empty. Then the following expansion of 
the minor of the matrix I + A can be obtained: 

(6.4.1) det (I + A) (l, ... , di, . . . , d r , ...,N\1, ...J 1 ,...J r , N) = 



(-)*=i 2^det A(/3 1 ,...,f3 r ,X\a 1 ,...,a r ,X) 

x 

Here the sum is over all subsets X of the set [1, . . . ,N], which do not intersect a and (5. 
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6.5. The minors of the matrix K. Applying two previous expressions it is easy to relate the 
minors of the matrix K defined in terms of the L matrix of the L-ensemble as 



(6.5.1) 



K = 1 - 



I + L 



with minors of the matrix L. Namely we have 
(6.5.2) 



, -r,, ,m det L(a, X\f3,X) 

where a = (cci, . . . , a r ), (3 = (/3i, . . . , (3 r ), a\ / /?_,-. Here again the sum is over all subsets X of the 
set [1, . . . , N]. 

6.6. Definition of PfafBan. The Pfaffian of a 2N x 2N antisymmetric matrix A = ||^4jfc||^ =1 is 
defined as 



(6.6.1) 



PiA= sgn(cr) A ili2 . . . A i2N _ ± 



*2JV 



(T=(u,...,i2Jv)eS2JV 
il<i2,...,l2N-l<i2N 
U<«3<---<*2JV-1 



6.7. PfafBan of the block matrix. The following formula is valid 

|A|(|A|-1) 



(6.7.1) 



Pf 



A 

-At 



det A 



6.8. The expansion of PfafBans of submatrices of the matrix J + A. Let A be a 2N x 2N 

antisymmetric matrix, 



(6.8.1) 



A[l',l",...,N',N"\l',l",...,N',N"} = 

A(l',l") ... A(1',N') A(1',N") 

-A(l',l") A(1",N') A(1",N") 



-A(1',N') -A(1',N") 
-A(1',N") -A(l",N") 







-A(N',N") 

and J be the antisymmetric matrix of format 2N x 2N defined by 



A(N',N") 




(6.8.2) 



J 



10 
-10 
1 
0-10 














1 
-1 



then the following expansion is valid (2 < 2k < N): 



Pf (J + A) 



l',2',...,2k'\l',2',...,2k' 



(6.8.3) 



X—(x' 1 ,x",...,x' d ,x'J l ) 
(x 1 ,...x d )c(2k+l,...,N) 



l",2",...,2k",X\l",2",...,2k",X 
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6.9. The Pfaffian of a submatrix of the inverse matrix. With A defined by equation (|6.8.1|) 
consider the inverse matrix B = A -1 . For 2 < 2m < N we find 



Pf B 



l / ,2 / ,...,2m'|l',2 / ,...,2m / 



1 



(6.9.1) x Pi A 



Pi A 

I", 2m", 2m + 1', 2m + l", . . . , N' , N"\ 

I", 2m", 2m + l', 2m + l", ...,N', N" 
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X 





X— \ Xq Xq 









Fig 1. The decomposition of the discrete set X 
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— \ Xq Xq X 



Fig 2. Unbalanced particle-particle configurations. The configuration X consists of four negative 
particles and one positive particle. The set Xq is chosen such that | 3£q I = — l-^+l- 
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Fig 3. Balanced particle-hole configurations. Z_ consists of three negative particles, \Z—\ = 3. Z + 
consists of one positive particle, \Z+\ = 1, and two holes in Xo, \Z+\ = 2 
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Splitting of X 


j£ — 1 3^_|_ ) — \ 1 1 3^o | 3^_|- ) 

|X | = S 


3C — 1 ^ — \ 3£q | 3£q 1 1 3£ | ) 


Point Configurations 


X = (X-\X+) = (xL u x H \x + ) 


Z = (Z \z+) = (z \z{ u zi 1 ) 




xi = x_n (x_ \ x ), x n = x_ n x 




L-ensembles 












Relations 
between configurations 


x£ = z_ 

X" = X \ Z{ 

x + = z» 


z_ = x£ 

Z| = X \ X u 

zi 1 = x + 



Table 1. L-ensemble and L-ensemble 
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L-ensemble 


A 7v-ensemble 


L-ensemble 


A Ar + 5-ensemble 


Splitting of 


X = (X_|X+) 

(X \ X 1 1 T IT ^ 

— yX— \ X.q U Xo | X_|_ ) 

|X | =S; |£+| =N 


X= (X_|X+) 


X= (X_|X+) 
(x \ t It iit ^ 

|X+| = N + S 


X= (X_|X+) 


Point 
configurations 


X = (X_\X + ) 

x+ = x n x+ 


X A = (X_\X+\X + ) 


Z=(Z_\Z + ) 

z_ = zr\X- 
z+ = z n x+ 


Z* = (Z_\X+\Z + ) 


Number of points 


arbitrary < N 


N 


arbitrary < AT + S 1 


N + S 


Weight 


h 


f 


h 


f = f 



Table 2. L-ensemble and polynomial ensemble 
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